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I .  INTRODUCTION 


A.  GENERAL 

Ship  silencing  continues  to  be  a  major  design  requirement 
in  the  construction  and  operation  of  submarines  and  surface 
ships.  It  has  been  long  established  that  both  equipment  and 
personnel  are  adversely  affected  by  unwatned  noise  and  vibra¬ 
tion.  In  today's  environment  both  subsurface  and  surface, 
with  recent  advances  in  acoustic  devices  the  very  surviva¬ 
bility  of  the  platform  is  directly  related  to  own  ship's 
noise.  Weapon  platforms  must  be  quiet  enough  to  escape  detec¬ 
tion  by  sophisticated  passive  sonar  devices  as  well  as  quiet 
enough  to  prevent  own  ship's  noise  from  interfering  with 
detection  and  prosecution  of  enemy  targets. 

The  more  traditional  approach  to  ship  silencing  has  been, 
and  continues  to  be,  that  of  vibrations  isolation.  This  ap¬ 
proach  requires  that  all  sources  of  vibration  to  be  placed 
on  resilient  energy-absorbing  mounts.  This  results  in  a 
significant  reduction  in  the  transmission  of  vibrational 
energy  from  the  equipment  and  main  engines  to  the  surrounding 
environment  but  is  not  totally  satisfactory  because  these  same 
mounts  must  serve  as  shock  mounts.  In  addition,  any  resilient 
mount  will  have  a  peak  efficiency  in  particular  frequency 
range  with  decreasing  efficiency  both  above  and  below  that 
frequency  range.  The  optimal  solution  would  be  to  use 
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resilient  mounts  in  conjunction  with  other  methods  of 
vibrational  control  to  reduce  the  amount  of  vibration  and 
noise  generated  as  well  as  reducing  the  amount  transmitted. 

As  with  any  design,  today's  ships  are  a  compromise  of 
requirements  and  cost  considerations.  However,  the  trend 
recently  has  been  toward  the  design  of  equipment  that  is 
inherently  quieter,  which  results  in  less  noise  and  vibra¬ 
tional  energy  at  the  source.  This  current  trend  in  design 
is  currently  being  approached  from  many  directions.  Among 
these  considerations  are: 

1.  Extreme  care  in  balancing  of  rotating  machinery  with 
consideration  of  vibrations  during  start  up  and  coast 
down  as  well  as  steady  state  operation. 

2.  Selection  of  operating  frequency  as  far  removed  from 
resonant  frequencies  as  possible. 

3.  Close  attention  to  component  tolerances. 

4.  Possible  use  of  high  damping  materials  for  machine 
elements  as  well  as  machine  casings  and  load  bearing 
structures . 

The  last  category,  use  of  high  damping  materials,  is  the 
furthest  from  being  utilized  to  its  maximum  potential  for 
noise  and  vibrational  reductions. 

Only  recently  has  the  designer  even  considered  the  inter¬ 
nal  capacity  of  a  material  alone  with  the  traditional  material 
properties  of  strength,  fatigue  resistance,  toughness  and 
corrosion  resistance. 
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There  are  several  commercially  available  high  damping 
alloys  that  meet  the  strength  requirements  for  most  shipboard 
applications.  It  is  recognized  that  each  of  the  candidate 
materials  may  have  unique  properties  that  may  cast  doubt  on 
their  actual  usefulness  in  the  United  States  Navy.  A  more 
important  problem  exists,  however.  That  is,  the  lack  of  a 
consistent  method  of  measuring  damping  for  a  material  under 
low  stress,  and  high  frequency  range. 

B.  DAMPING 

"Structural  damping"  refers  to  a  structure's  or  struc¬ 
tural  component’s  capacity  for  dissipating  energy,  or,  more 
precisely,  to  its  capacity  for  removing  from  a  structural 
vibration  some  of  the  energy  associated  with  that  vibration. 
The  energy  removed  may  be  converted  directly  into  heat,  trans¬ 
ferred  to  connected  structures  or  ambient  media  [Ref.  1] . 

Damping  has  two  primary  effects:  (1)  It  limits  the 
steady-state  motions  of  structures  or  systems  in  situations 
where  these  motions  are  controlled  by  an  energy  balance;  and 
(2)  It  increases  the  rates  at  which  the  free  (i.e.,  unforced) 
vibrations  of  structures  decay. 

Consider,  for  example,  the  classical  lumped-parameter 
mass-spring  dashpot  system  driven  by  a  steady  sinusoidal  force 
that  acts  on  the  mass.  For  such  a  system  it  is  possible  to 
make  the  following  observations:  (1)  For  excitation  frequen¬ 
cies  that  are  considerably  lower  than  the  system  resonance 


frequency  the  applied  force  essentially  is  balanced  by  the 
spring  force.  The  mass  and  dashpot  here  have  virtually  no 
effect.  (2)  For  excitation  frequencies  that  are  consider¬ 
ably  higher  than  the  system  resonance  frequency,  the  applied 
force  essentially  is  opposed  by  only  the  inertia  force  of  the 
mass.  In  this  case  the  spring  and  dashpot  have  virtually  no 
effect. 

(3)  For  excitation  frequency  near  the  system  resonance 
frequency,  the  spring  force  and  the  inertia  force  essentially 
cancel  each  other,  leaving  only  the  dashpot  (i.e.,  the  damping 
element)  to  oppose  the  externally  applied  force.  One  may 
note  that  changes  in  the  damping  do  not  affect  the  system 
response  for  the  first  two  of  the  above  cases,  but  for  an 
oscillatory  force  acting  at  the  system's  resonance,  the  vibra¬ 
tion  amplitude  decreases  with  increasing  damping. 

In  addition  to  the  case  of  resonant  or  near- resonant 
excitation,  there  occur  several  other  important  steady-state 
or  near-steady-state  situations  in  which  the  responses  of 
systems  or  structures  are  controlled  by  a  balance  between  the 
energy  input  and  the  energy  dissipations.  These  situations 
include  cases  of  (1)  broad-band  excitation,  and  (2)  spatially 
periodic  excitation  where  the  spatial  period  matches  that  of 
freely  propagating  waves.  Because  broad-band  excitation  gener¬ 
ally  encompasses  several  structural  resonance  frequencies, 
the  structural  responses  of  the  excited  modes,  and  each 
resonant  modal  response  is  controlled  by  damping. 
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If  a  structure  (or  a  mechanical  system)  is  deflected  from 


its  equilibrium  position  and  then  released,  the  structure 
vibrates  with  ever-decreasing  amplitude  on  the  result  of 
damping,  i.e.,  as  the  result  of  energy  being  removed  from 
the  oscillatory  motions.  Greater  damping  corresponds  to  the 
dissipation  per  cycle  of  a  greater  fraction  of  the  vibratory 
energy,  thus  resulting  in  more  tapid  decay  of  the  vibrations. 

In  a  somewhat  similar  manner,  increased  damping  also  results 
in  the  more  rapid  decay  of  freely  traveling  waves.  For  exam¬ 
ple,  if  a  long  beam  is  subjected  to  an  oscillatory  transverse 
force  of  a  constant  amplitude  and  frequency,  such  as  near  the 
center  of  the  beam,  then  flexural  waves  travel  away  from  the 
driving  point  in  both  directions.  As  the  result  of  damping, 
the  amplitude  of  these  waves  decreases  with  increasing  distance 
from  the  point  of  application — with  greater  damping  leading 
to  lesser  amplitudes  at  a  given  distance. 

The  practical  consequences  of  effects  mentioned  above 
generally  are  the  reasons  one  is  interested  in  systems  or 
devices  that  increase  damping.  Reductions  in  resonant  or 
random  responses  result  in  decreased  oscillatory  stresses 
and  increase  in  the  fatigue  life  of  structures,  in  the  relia¬ 
bility  of  mechanical  devices,  and 'in  mechanical  impedance 
(which  tends  to  improve  the  effectiveness  of  vibration 
isolation) . 

Reduction  of  the  spatially  resonant  responses  of  a  wall 
or  panel  leads  to  decreased  sound  transmission  through  that 
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structure  for  frequencies  above  the  coincidence  frequency. 
Increased  attenuation  of  propagating  waves  results  in  lesser 
transmissions  of  vibrations  to  neighboring  structures.  More 
rapid  decay  of  free  vibrations  reduces  the  "ringing"  sound  of 
structures,  thus  leading  to  less  noise,  particularly  from 
structures  excited  by  repetitive  impacts .  This  also  tends 
to  reduce  structural  fatigue. 


C.  MEASURES  OF  DAMPING 

Damping  may  be  quantified  in  terms  of  any  of  the  previously 
discussed  primary  effects .  The  corresponding  commonly  used 
measures  of  damping,  defined  below,  are  interrelated  as 
follows:  [Refs.  2,3] 


4>  _  2.20  _  At  _  5  _  Ax 

27  "  T~,n  *  2T.'3  r  =  7  *  T777 

nod  n 


(1) 


The  loss  factor  n  and  damping  capacity  ^  are  defined 
directly  in  terms  of  the  cyclic  energy  dissipation;  the 
damping  capacity  represents  the  fraction  of  the  system's 
vibrational  energy  that  is  dissipated  per  cycle  of  the  vibra¬ 
tion,  and  the  loss  factor  similarly  is  defined  as  the  frac¬ 
tion  of  the  system's  energy  that  is  dissipated  per  radian  of 
the  vibratory  motion. 

On  the  other  hand,  Tgg,  are  related  to  the  rate  of 
decay  of  free  vibrations.  TgQ  denotes  the  reverberation 
time  (in  seconds) ,  defined  (in  analogy  to  the  related 
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room-acoustic  measure)  as  the  time  within  which  the  vibra¬ 
tion  level  of  a  system  vibrating  freely  at  a  frequency 
fn  (Hz)  decreases  by  60  dB  (i.e.,  the  amplitude  decreases  to 
1/1000  of  its  original  value) .  A  related  measure,  decay  rate 
(dB/sec) ,  represents  the  rate  of  reduction  of  the  vibra¬ 
tion  (acceleration  or  displacement  velocity)  level.  The 
logarithmic  decrement  6  is  defined  as  the  natural  logarithm 
of  the  ratio  of  a  peak  excursion  of  a  freely  vibrating  system 
to  the  peak  excursion  one  cycle  (period)  later.  The  spatial 
decay  A  (dB/wavelength)  represents  the  reduction  in  the 
steady-state  vibration  level  with  distance  that  occurs  along 
a  long  beam  vibrating  in  flexure. 

It  should  be  noted  that  none  of  the  measurements  of  damping 
depend  on  how  the  energy  is  dissipated.  Within  a  cycle  these 
measures  make  no  reference  to  any  damping  mechanism.  On  the 
other  hand,  some  other  commonly  employed  measures  of  damping 
are  defined  on  the  basis  of  viscous  damping  (i.e.,  damping 
that  results  from  a  retarding  force  that  is  proportional  to 
the  velocity.  The  ratio  of  the  magnitude  of  that  force  to 
the  velocity  is  called  the  viscous  damping  coefficient  and 
is  commonly  designated  by  c. 

If  a  simple  mass-spring-dashpot  system  (where  the  dashpot 
provides  a  viscous  retarding  force  characterized  by  c)  is 
deflected  from  equilibrium  and  released,  it  typically  oscil¬ 
lates  with  ever  decreasing  amplitude.  However,  if  c  is  made 
large  enough,  no  oscillations  occur.  Instead,  the  system 


creeps  toward  its  equilibrium  position,  never  traversing  it. 

The  smallest  viscous  damping  coefficient  for  which  this  non- 
oscillatory  behavior  is  obtained — i.e.,  the  viscous  damping 
that  represents  the  dividing  line  between  oscillatory  and 
nonoscillatory  behavior  is  called  the  critical  (viscous) 
damping  coefficient  c  .  It  obeys  c  =  2 /km  where  m  denotes 
the  mass  and  k  the  spring  stiffness.  The  "damping  ratio" 
c/c  ,  also  called  the  fraction  of  critical  deunping  and  often 

G 

given  in  terms  of  "percent  of  critical  damping",  is  widely 
used  to  indicate  deunping  magnitudes. 

Two  other  measures  of  deunping  are  derived  from  the  steady 
state  behavior  of  an  ideal  linear  mass-spring-dashpot  system 
that  is  driven  by  a  sinusoidal  force  of  constant  amplitude 
[Ref.  4] .  The  amplification  at  resonance,  often  called  "the 
Q"  of  the  system  is  defined  as  the  ratio  of  the  amplitude  that 
results  at  resonance  to  the  amplitude  that  is  obtained  if 
the  force  acts  quasi-statically  (i.e.,  a  frequency  consider¬ 
ably  below  resonance) .  The  proportional  bandwidth  b  takes 
account  of  the  damping  related  broadening  of  the  peak  in  a 
plot  of  response  amplitude  versus  frequency;  this  nondimens ional 
bandwidth  is  defined  as  f/fR,  where  f  denotes  the  difference 
between  the  two  frequences  (one  above,  and  one  below  the 
resonance  frequency  fn)  at  which  the  square  of  the  response 
amplitude  is  one-half  of  its  maximum  value.  For  values  of 
damping  below  critical,  the  aforementioned  measure  of  damping 
are  related  to  each  other  and  to  the  previously  discussed 
loss  factor  as: 
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D.  DAMPING  MECHANISMS 

For  linear,  viscously  damped  systems,  all  of  the  measures 
of  damping  discussed  previously  are  independent  of  the  ampli¬ 
tude.  Amplitude  independence  also  occurs  for  other  deunping 
mechanisms,  and  approximate  amplitude-independence  occurs  for 
almost  all  systems,  provided  the  damping  is  relatively  small. 
Thus,  amplitude  independence  cannot  be  taken  as  an  indication 
that  a  system  is  viscously  damped;  nevertheless,  damping  of 
systems  with  unknown  energy  dissipation  mechanism  is  often 
characterized  in  terms  of  an  equivalent  viscous  damping 
coefficient. 

Indeed,  much  analysis  is  carried  out  with  the  (usually 
tacit)  assumption  that  the  damping  is  viscous — largely  be¬ 
cause  viscous  damping  leads  to  linear  differential  equations 
that  can  be  solved  relatively  easily.  It  is  fortunate  that 
for  many  practical  problems,  e.g.,  where  only  certain  response 
maxima  are  of  concern,  the  details  of  the  damping  mechanism 
are  unimportant.  Thus,  one  may  obtain  reasonable  response 
predictions  even  with  inaccurate  damping  force-versus-velocity 
representations.  However,  realistic  damping  models  are  re¬ 
quired  for  the  analysis  of  cases  where  one  is  interested  in 
details  of  the  response  motions. 


Unlike  mass,  a  single  physical  phenomenon  and  stiffness, 
which  results  from  a  very  few  physical  effects,  damping  may 
be  caused  by  a  great  variety  of  phenomena.  These  phenomena 
include  mechanical  hysteresis  (also  called  material  damping 
or  internal  friction) ,  electromagnetic  effects  (notably 
eddy  currents) ,  friction  due  to  motion  relative  to  fluids 
or  solid  surfaces,  and  energy  transport  to  adjacent  structural 
components  or  fluids  (including  by  acoustic  radiation) . 

This  great  variety  of  phenomena  that  can  produce  damping 
generally  make  damping  difficult  to  predict  and  to  eliminate, 
but  enables  one  to  conceive  a  variety  of  means  for  increasing 
it. 
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II.  NATURE  OF  THE  PROBLEM 


A.  BACKGROUND 

To  date,  the  majority  of  damping  research  has  been  con¬ 
ducted  on  test  specimens  that  are  subjected  to  high  stresses 
(i.e.,  torsional  vibrations  and  vibrating  cantilevered  beams). 
For  most  United  States  Naval  applications,  structures  and 
components  are  designed  for  operation  at  relatively  low 
stresses. 

Damping  characteristics  are  shown  to  be  highly  stress 
dependent  [Ref.  5].  Therefore,  the  damping  values  obtained 
at  high  stresses  do  not  generally  apply  for  most  naval 
applications  of  ship  silencing  problems. 

David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(NSRDC) ,  Annapolis,  Maryland  is  the  focal  point  for  testing 
and  evaluating  new  high  damping  material  for  possible  use 
for  the  United  States  Navy.  Realizing  the  wide  scope  of  the 
work  necessary  for  the  testing  and  evaluation  of  any  new 
material,  the  author  decided  to  approach  NSRDC  with  a  proposal 
to  offer  assistance  in  determining  the  damping  of  candidate 
material. 

B.  OBJECTIVE 

To  design  a  test  procedure  that  allows  the  measurement  of 
damping  in  a  plate  specimen  at  low  stress  with  the  following 
variables: 
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1.  Plate  specimen  (40  inches  x 14  inches  x l  inch)  and 

40  inches  x 14  inches  x  2  inches) ,  or  in  proportionally 
reduced  sizes. 

2.  Frequency:  Damping  testing  will  be  conducted  in  the 
acoustic  range  (100-20,000  Hz). 

3.  Environment:  Testing  will  be  conducted  in  lab  air 
and  in  nondistilled  water.  Temperature  will  range 
from  30°F  to  90°F. 

The  test  procedure  includes  a  complete  modal  analysis 
that  is  expandable  to  other  geometric  shapes.  Actual  testing 
and  verification  of  the  procedure  are  conducted  on  the 
Hewlett-Packard  5451c  Digital  Fourier  Analysis  System. 

General  information  on  the  HP-54 51C  system  is  included  in 
Appendix  A. 

C.  SCOPE  OF  WORK  TO  BE  COMPLETED 

NSRDC  has  provided  four  specimens  for  testing  and  to 
establish  baseline  data  for  future  work.  The  test  specimens 
currently  on  hand  are:  A)  Cast  manganese  bronze,  code  DEQ; 

B)  Cast  nickel-aluminum  bronze,  code  FTC;  C)  Steel  plate, 
HY-130,  code  FTW;  and  D)  Aluminum  alloy,  5086-H116,  code  ESX. 
It  is  anticipated  that  after  baseline  data  has  been  compiled 
for  both  air  and  water  tests  that  additional  specimens  will 
be  provided  by  NSRDC.  These  additional  specimens  are  ex¬ 
pected  to  include  high  damping  alloys  such  as  Sonoston  and 
Incramute,  as  well  as  constrained  layer  specimens  and  com¬ 
posites,  both  metallic  and  non-metallic. 
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For  the  purpose  of  this  thesis  the  following  work  was 
completed:  (A)  design  and  construction  of  the  test  chamber, 

(B)  design  of  a  test  procedure  for  the  evaluation  of  damping 
of  specimen  provided  by  NSRDC,  (C)  complete  tank  character¬ 
ization  by  impulse  hammer  techniques,  (D)  measurement  of 
damping  of  the  cast  nickel-aluminum  bronze,  code  FTC,  speci¬ 
men  by  impulse  hammer  technique.  All  of  the  above  testing 
was  completed  in  lab  air  at  normal  ambient  temperature. 
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III.  DESIGN  OF  TEST  CHAMBER  AND  THEORY  OF  CHARACTERIZATION 


To  accomplish  the  aforementioned  objectives  four  major 
steps  must  be  accomplished:  (A)  design  of  a  test  chamber, 

(B)  characterization  of  the  test  chamber,  (C)  design  of 
experimental  procedure  of  the  measurement  of  damping,  and 
(D)  utilization  of  the  above  test  procedure  for  damping 
measurement  on  a  plate  specimen  supplied  by  NSRDC. 

A.  DESIGN  OF  TEST  CHAMBER 

Design  of  the  test  chamber  was  based  on  three  major 
considerations;  (1)  size  of  the  largest  specimen,  (2)  weight 
of  the  largest  specimen,  and  (3)  environmental  conditions 
desired  for  testing.  In  addition,  the  fixture  to  support 
the  specimen  had  to  hold  the  specimen  rigid  so  that  there 
would  be  no  swinging. 

The  overall  outside  dimension  of  the  chamber  (Figure  1) 
was  chosen  to  be  a  rectangular  box  48  inches  wide  by  36  inches 
deep  by  72  inches  high.  This  allowed  for  ample  clearance  for 
the  largest  specimen  (40  inches  *  14  inches  *  2  inches)  on  all 
surfaces  and  to  limit  the  amount  of  reflected  acoustic  energy 
if  acoustic  absorption  material  was  required  on  the  inside 
surfaces.  The  test  chamber  was  constructed  from  one  quarter 
inch  plate  steel  supported  by  a  frame  constructed  from  2  inch 
steel  thick  wall  square  tubing.  The  2  inch  square  tubing 


Figure  1.  Experiment  test  chamber 


was  welded  at  all  mating  edges  to  give  the  desired  rigidity 
and  to  prevent  water  entrapment  (future  corrosion  problems) . 

The  specimen  fixture  (Figure  2)  was  constructed  from 
3  inches  x  2  inches  thick  wall  steel  tubing  and  is  designed 
to  give  a  deflection,  from  a  400  pound  specimen,  of  less  than 
0.04  inch. 

To  accommodate  various  sized  specimen  the  random  frequency 
exciter  base  (Figure  3)  can  be  adjusted  within  the  2  inch 
rectangular  tubing  frame  work  of  the  upper  half  of  the  rear 
panel  (Figure  4) .  The  random  acoustic  exciter  can  be  adjusted 
from  8  inches  below  the  specimen  fixture  to  24  inches  below 
the  specimen  fixture.  This  feature  allows  for  a  wide  variety 
of  centered  excitation  from  the  end  of  a  specimen  plate  to 
the  center  of  a  specimen  plate. 

To  allow  the  removal  and  placement  of  various  size  speci¬ 
men,  the  front  panel  of  the  test  chamber  is  bolted  to  the 
test  chamber  frame  and  a  water  gasket  is  formed  by  using  a 
RTV  silicone  adhesive  sealant.  Because  of  the  bonding  of 
the  RTV  all  water  test  must  be  conducted  with  the  steel  front 
panel.  For  air  testing  the  optional  clear  front  panel  may 
be  used  because  no  sealant  is  required. 

The  preservation  of  the  test  chamber  was  accomplished  by 
sandblasting  all  steel  parts,  after  protecting  front  panel 
and  specimen  fixture  bolt  threads,  one  spray  paint  coat  of 
zinc  chromate  primer  and  one  hand  brushed  top  coat  of  enamel. 

To  allow  testing  in  lab  air  and  nondistilled  water,  the 
test  chamber  is  water  tight  with  a  1  inch,  valved,  drain 


Ficjuro  4.  Inside  rear  panel  of  test  chamber 


line  near  the  bottom  of  the  tank.  Heating  and  cooling  of  the 
tank  is  accomplished  by  copper  tubing  (3/8  inch  O.D.)  coiled 
over  entire  outside  surface  of  the  back  panel  (Figure  5) . 

The  coils  are  spaced  3  inches  apart  and  are  40  inches  wide. 

To  heat  or  cool  the  test  chamber  a  mixture  of  50%  ethol  glycol 
(automative)  antifreeze  is  supplied  from  a  storage  tank 
(Figure  6)  at  the  desired  temperature  and  is  pumped  through 
the  heating  and  cooling  coil. 

The  entire  test  chamber  is  insulated  with  one  and  one 
half  inches  of  rigid  foam  insulation,  with  the  exception  of 
the  rear  panel,  which  is  insulated  with  standard  residential 
batt  insulation. 

B.  THEORY  OF  THE  CHARACTERIZATION  OF  TEST  CHAMBER 

The  characterization  of  the  test  chamber  was  accomplished 
by  use  of  the  impulse  technique  for  structural  frequency 
response  testing  [Ref.  7]  and  the  Fourier  transfer  function 
capability  of  the  HP-54 51C  Fourier  Analyzer. 

1 .  Theory  of  Frequency  Response  Function 

The  measurement  of  the  frequency  response  function  is 
the  heart  of  modal  analysis.  The  frequency  response  function 
H (f )  is  defined  in  terms  of  the  single  input/single  output 
system,  as  the  ratio  of  the  Fourier  transforms  of  the  system 
output  or  response  v(t)  to  the  system  input  or  excitation 
u(t) ,  Equation  (3). 


H  ( f ) 


V(f) 

uTFT 


(3) 
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Figure  6.  Cooling/heating  water  tank  with  pump,  heater  and  controller 


where : 


V(f )  =  Fourier  transform  of  system  output  v‘(t) 

U(f)  =  Fourier  transform  of  system  input  u(t). 

The  only  requirements  for  a  complete  description  of 
the  frequency  response  function  are  that  the  input  and  output 
signals  be  Fourier  transformable,  a  condition  that  is  met  by 
all  physically  realizable  systems,  and  that  the  input  signal 
be  non-zero  at  all  frequencies  of  interest.  If  the  system 
is  nonlinear  or  time-variant,  the  frequency  response  function 
will  not  be  unique,  but  will  be  a  function  of  the  amplitude 
of  the  input  signal  in  the  case  of  a  nonlinear  system  and 
a  function  of  time  in  the  case  of  a  system  with  time-varying 
properties. 

The  frequency  fesponse  function  may  be  computed  directly 
from  the  definition  as  the  ratio  of  the  Fourier  transforms 
of  the  output  and  input  signals.  However,  better  results  are 
obtained  in  practice  by  computing  the  frequency  response 
function  as  the  ratio  of  the  cross-spectrum  between  the  input 
and  output  to  the  power  spectrum  of  the  input.  Equation  (4). 

This  relationship  is  derived  by  multiplying  the  numerator  and 
denominator  of  the  right-hand  side  of  Equation  (1)  by  the 
complex  conjugate  of  the  input  Fourier  transform. 

Guv(f) 

H(f)  =  (4) 
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where : 


Guv(f) 

=  U*(f)V(f), 

cross-spectrum 

between  u(t) 

and  v(t) 

Guv(f> 

=  U* (f ) U(f ) , 

power  spectrum 

of  u(t) 

u* 

=  complex  conjugate  of  U(f). 

The  usefulness  of  this  form  of  the  frequency  response  function 
can  be  seen  by  considering  the  practical  single  input/single 
output  measurement  situation,  where  m(t)  and  n(t)  represent 
noise  at  the  input  and  output  measurement  points,  respectively. 

The  measured  frequency  response  function  H'(f)  is 
given  by  the  equation: 


H'  (f) 


Y(f)  _  V(f)  +  N(f) 

xTfT  Off)  +  MtfJ 


(5) 


where  the  upper  case  letters  denote  the  Fourier  transform  of 
the  corresponding  time  domain  signals. 

In  this  form,  the  measured  frequency  response  will  be 
a  good  approximation  of  the  true  frequency  response  only  if 
the  measurement  noise  at  both  the  input  and  output  measure¬ 
ment  points  is  small  relative  to  the  input  and  output 
signals .  Multiplying  the  numerator  and  denominator  of  the 
right-hand  side  of  Equation  (5)  by  the  complex  conjugate  of 
X(f)  yields 


H *  (f ) 


Guv(f)  *  Gun(f)  +  Gmv(f)  +  Gmn(f> 
6u(fT  +  Gu»(f>  +  Gmn<f)  +  Gm(f) 


(6) 
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Now,  if  the  measurement  noise  signals  m(t)  and  n(t) 
are  noncoherent  with  each  other  and  with  the  input  signal 
u(t) ,  then  the  expected  value  of  the  cross-spectrum  terms 
involving  m  and  n  in  Equation  (6)  will  equal  zero,  yielding 


H  *  ( f ) 


Guv<f) 

Gu(f)  +Gm(f) 


1  + 


H(f) 

G  ( f ) 
m 

G  (f) 
u 


(7) 


where  H(f)  is  the  desired  true  frequency  response  function. 

Thus,  if  the  noise-to-signal  ratio  at  the  input 
measurement  point  [Gm( f )/Gu(f ) ]  is  much  less  than  1,  the 
measured  frequency  response  will  closely  approximate  the 
desired  true  frequency  response  function. 

It  should  be  pointed  out  here  that  there  is  an  inher¬ 
ent  bias  error  associated  with  the  computation  of  the  cross¬ 
spectrum  and  the  magnitude  of  this  bias  error  is  inversely 
proportional  to  the  number  of  averages  in  the  computation. 

Thus,  the  greater  the  measurement  noise,  the  greater  the  num¬ 
ber  of  averages  required  to  approach  the  expected  value  of 
the  cross- spectrum  between  the  input  and  the  output  measure¬ 
ment  signals.  With  measurement  techniques  employing  many 
averages,  the  bias  error  can  usually  be  reduced  to  an  insigni¬ 
ficant  level  so  that  it  is  only  necessary  to  minimize  the  noise 
in  the  measurement  of  the  input  signal.  However,  if  there  is 
significant  measurement  noise  and  only  a  few  averages  are 
used,  then  the  computed  values  of  the  cross -spectrum  terms 
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involving  the  noise  signals  in  Equation  (6)  can  be  large 
relative  to  the  true  cross-spectrum,  with  resulting  large 
errors  in  the  measured  frequency  response  function.  In 
general,  only  a  few  averages  are  used  in  the  impulse  technique; 
otherwise,  one  of  its  major  advantages — its  speed — is  lost. 
Therefore,  it  is  important  to  minimize  measurement  noise  in 
both  the  input  and  output  signals  when  using  the  impulse 
technique.  The  cross- spectrum  bias  error  and  its  effects 
are  discussed  in  more  detail  in  Reference  8 . 

Coherence  Function.  There  is  another  important  reason  for 
computing  the  frequency  response  function  in  terms  of  the 
cross- spectrum:  it  allows  the  computation  of  the  coherence 

function  between  the  input  and  output  signals.  The  coherence 
function  is  defined  by  the  equation 


2  lGxy(f)!2 

Yxyir>  Gx(f)Gy(f) 


(8) 


According  to  the  definitions  of  the  power  spectrum 
and  the  cross- spectrum,  the  coherence  function  will  be  iden¬ 
tically  equal  to  1  if  there  is  no  measurement  noise  and  the 
system  is  linear.  The  minimum  value  of  the  coherence  function, 
which  occurs  when  the  two  signals  are  totally  uncorrelated,  is 
0.  Thus,  the  coherence  function  is  a  measure  of  the  con¬ 
tamination  of  the  two  signals  in  terms  of  noise  and  nonlinear 
effects,  with  very  low  contamination  indicated  for  values 


close  to  1 


Since  the  cross-spectrum  is  included  in  the  definition 
of  the  coherence  function,  the  cross-spectrum  bias  error  must 
Le  reduced  to  an  acceptable  level  if  a  good  statistical  esti¬ 
mate  of  the  coherence  function  is  to  be  achieved.  As  stated 
above,  the  number  of  averages  used  in  the  impulse  technique 
is  usually  not  great  enough  to  significantly  reduce  the  bias 
error.  However,  the  coherence  function  is  still  useful  for 
indicating  the  importance  of  noise  in  the  impulse  technique. 
This  is  because  noise  in  the  signals  causes  variance  in  the 
value  of  the  coherence  function  with  frequency.  This  effect 
is  illustrated  in  the  section  on  measurement  procedures. 

2 .  Display  of  Frequency  Response 

The  frequency  response  function  is  complex — that  is, 
it  has  associated  with  it  both  magnitude  and  phase.  There¬ 
fore,  it  can  be  displayed  in  a  number  of  forms,  including 
magnitude  and  phase  versus  frequency,  real  and  imaginary 
magnitudes  versus  frequency,  and  imaginary  magnitude  versus 
real  magnitude.  Each  of  these  types  of  displays  has  its 
own  particular  usefulness.  The  most  common  type  of  display 
for  structural  frequency  response  data  is  magnitude  and  fre¬ 
quency  plotted  logarithmically.  This  type  of  display,  with 
the  magnitude  in  terms  of  compliance  (ratio  of  displacement 
to  force)  is  called  a  Bode  plot.  In  this  form  of  the  fre¬ 
quency  response  function,  resonances  occur  as  peaks  in  com¬ 
pliance  plots  (points  of  maximum  dynamic  weakness)  and  all 
resonance  peaks  of  equal  damping  have  the  same  width  regardless 
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of  resonance  frequency.  Lines  of  constant  dynamic  stiffness 
have  zero  slope,  and  mass-dominated  frequency  response  lines 
have  a  -12  dB-per-octave  slope.  Resonances  occur  as  nearly 
circular  arcs  in  the  complex  plane  (real  versus  imaginary 
plot)  with  frequency  increasing  in  a  clockwise  direction  around 
the  arc.  In  the  case  of  real  normal  modes  (which  occur  in 
systems  with  relatively  low  damping  and  with  resonances  well- 
separated  in  frequency) ,  each  resonance  arc  is  approximately 
tangent  with,  and  lies  below,  the  real  axis  and  is  symmetric 
about  the  imaginary  axis  when  the  frequency  response  is  ex¬ 
pressed  as  compliance.  The  complex  plane  plot  is  useful  when 
certain  types  of  analytical  curve  fitting  operations  are 
being  performed  on  the  frequency  response  data.  The  plots 
of  the  real  and  imaginary  magnitudes  of  frequency  response 
versus  frequency  are  most  useful  when  dealing  with  real  normal 
modes.  In  this  case  the  resonances  will  occur  as  peaks  in 
the  imaginary  magnitude  plot  and  the  real  magnitude  will  pass 
through  zero  at  the  resonance  frequency  when  the  frequency 
response  is  expressed  as  compliance. 

The  frequency  response  characteristics  of  a  structural 
element  are  determined  by  measuring  a  set  of  cross-frequency 
response  functions  as  discussed  in  Reference  9.  The  cross¬ 
frequency  response  functions  may  be  obtained  by  exciting  at 
one  location  on  the  structure  and  measuring  response  at 
various  locations,  or  by  measuring  the  response  at  a  single 
location  to  excitation  at  various  locations.  The  resulting 


frequency  response  functions  comprise  one  column  of  the  trans¬ 
fer  matrix  in  the  first  case,  and  one  row  of  the  transfer 
matrix  in  the  second  case.  Either  set  will,  in  general,  com¬ 
pletely  define  the  modal  characteristics  of  the  structural 
element.  In  mathematical  terms  the  set  of  frequency  response 
functions  yields  the  eigenvalues  and  eigenvectors,  which  are, 
in  general,  complex  terms.  The  real  part  of  an  eigenvalue  is 
the  damping  and  the  imaginary  part  is  the  frequency  associated 
with  a  given  resonance.  Each  eigenvector  defines  a  resonance 
mode  shape. 

With  real  normal  modes,  each  point  on  a  structure  is 
either  exactly  in-phase  or  exactly  180  degrees  out-of-phase 
with  any  other  point  at  the  resonance  frequency.  Certain 
types  of  damping  which  are  often  encountered  in  practice  will 
cause  the  eigenvectors  to  have  non-zero  imaginary  components, 
resulting  in  complex  mode  shapes.  When  a  mode  is  complex, 
the  relative  phase  associated  with  a  point  on  a  structure  is 
some  value  other  than  0  or  180  degrees,  with  the  result  that 
node  lines  (lines  of  zero  deflection)  are  not  stationary. 
Precise  description  of  complex  modes  requires  that  some  type 
of  analytical  curve  fitting  technique  be  applied  to  the  fre¬ 
quency  response  data. 

The  frequency  response  function  of  an  operating  system 
can  be  computed  if  the  system  input  and  output  signals  meet 
previously  stated  requirements  of  Fourier  transformability 
and  non-zero  value,  assuming  the  system  input  and  response 
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can  be  measured.  However,  in  practice  there  are  usually 
multiple  inputs  to  the  system — either  several  inputs  at 
different  locations  or  inputs  in  more  than  one  direction  at 
a  given  location.  In  the  case  of  multiple  coherent  inputs, 
the  complexity  of  the  analysis  is  greatly  increased.  For 
this  reason,  and  the  difficulty  of  accurately  monitoring 
operating  inputs,  frequency  response  measurements  are  usually 
made  by  applying  the  system  input  "artificially"  through 
some  type  of  exciter.  It  is  in  the  form  of  the  input  signal 
and  the  way  it  is  applied  to  the  structure  that  the  wide 
variety  of  frequency  response  testing  techniques  arises. 

The  usefulness  of  the  impulse  technique  lies  in  the 
fact  that  the  energy  in  an  impulse  is  distributed  continu¬ 
ously  in  the  frequency  domain  rather  than  occurring  at  dis¬ 
crete  spectral  lines  as  in  the  case  of  periodic  signals. 

Thus,  an  impulse  force  will  excite  all  resonances  within  its 
useful  frequency  range.  The  extent  of  the  useful  frequency 
range  of  an  impulse  is  a  function  of  the  shape  of  the  im¬ 
pulse  and  its  time  duration.  For  a  square  pulse  the  frequen¬ 
cies  of  the  zero  crossings  are  at  integral  multiples  of  the 
inverse  of  the  time  duration  of  the  impulse,  illustrating 
the  very  important  inverse  relationship  between  the  time 
duration  of  an  impulse  and  its  frequency  content. 

The  useful  frequency  range  of  an  impulse  is  also  a 
function  of  the  shape  of  the  impulse.  By  varying  the  weight 
and  hardness  of  an  impacting  device  and  the  manner  in  which 
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the  impact  is  applied,  the  shape  and  time  duration  of  the 
impulse  produced  can  be  varied  to  suit  the  measurement 
requirements . 

a.  Nonlinearities  in  the  Structures 

Excitation  of  a  nonlinear  system  by  a  pure-random 
signal  will  yield  the  best  estimate  (in  a  mean-square  sense) 
of  the  linear  system  response.  Excitation  by  a  pure  sine 
wave  is  also  useful  for  studying  nonlinear  systems  because 
it  allows  precise  control  of  the  input  spectrum  level.  How¬ 
ever,  the  impulse  technique,  because  of  its  very  high  ratio 
of  peak  level  to  total  energy,  is  particularly  ill-suited 
for  testing  nonlinear  systems.  Therefore,  it  is  important 
to  understand  the  various  types  of  nonlinearities  that  can 
occur  in  structural  systems  and  to  be  able  to  recognize 
nonlinearities  in  measured  frequency  response  functions. 

One  of  the  most  common  types  of  nonlinearities 
encountered  in  structures  is  that  due  to  clearance  between 
parts.  This  type  of  nonlinearity  is  frequently  encountered, 
for  example,  when  testing  gear  systems  and  shafts  mounted  in 
bearings.  The  effects  of  this  type  of  nonlinearity  on  measured 
frequency  response  functions  when  using  impulse  excitation 
are  poor  estimates  of  static  stiffness  values  and  poor  re¬ 
peatability  of  the  frequency  response  estimates.  Also,  the 
apparent  damping  in  the  estimates  will  be  greater  than  the 
actual  examples. 
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The  best  method  of  dealing  with  this  type  of 
nonlinearity  is  to  preload  the  system  to  take  up  clearances. 
Care  must  be  taken  when  this  is  done,  however,  because  any 
preload  will  change  the  boundary  conditions  of  the  structure 
and  can  itself  lead  to  erroneous  frequency  response  esti¬ 
mates.  The  usual  approach  is  to  apply  the  preload  through  a 
very  soft  spring  so  that  the  resonances  associated  with  the 
preload  lie  below  the  frequency  range  of  interest. 

Another  type  of  nonlinearity  that  is  frequently 
encountered  is  nonlinear  damping.  Nonlinear  damping  effects 
are  usually  associated  with  joints  in  the  structure,  where 
the  damping  is  a  function  of  the  relative  displacement  at  the 
joint.  In  general,  the  frequency  response  estimates  obtained 
by  the  impulse  technique  will  agree  most  closely  with  those 
obtained  with  a  low  level  of  continuous  excitation.  However, 
if  the  point  of  excitation  is  close  to  a  location  where  non¬ 
linear  damping  occurs,  there  will  be  high  relative  motion  at 
that  location,  and  the  apparent  damping  in  the  measured  fre¬ 
quency  response  will  be  high.  In  systems  with  low  damping, 
this  will  give  the  measured  frequency  response  a  discontinuous 
appearance,  due  to  the  varying  level  of  damping  as  the 
response  to  the  impulse  attenuates  with  time. 

The  third  type  of  nonlinearity  that  commonly 
occurs  in  structures  is  load-sensitive  stiffness,  where  the 
spring  rate  of  elastic  elements  either  increases  or  decreases 
with  load.  The  most  direct  way  to  identify  this  type  of 


nonlinearity  is  to  measure  frequency  response  as  a  function 
of  static  preload  and  observe  the  change  in  resonance 
frequencies. 

b.  Signal  Processing 

The  particular  characteristics  of  an  impulsive 
force  signal  and  the  resulting  structural  response  signal 
make  the  impulse  technique  especially  susceptible  to  two 
problems:  noise  and  truncation  errors.  While  these  problems 

occur  to  some  extent  with  other  frequency  response  testing 
techniques,  their  unique  importance  in  the  impulse  technique 
requires  special  signal  processing  methods. 

It  was  pointed  out  in  the  previous  section  that 
the  usable  frequency  range  for  an  impulse  depends  on  the 
shape  and  time  duration  of  the  impulse.  In  order  to  insure 
that  there  is  sufficient  force  over  the  frequency  range  of 
interest,  it  is  necessary  that  the  first  zero  crossing  of 
the  Fourier  transform  of  the  impulse  be  well  above  the  maxi¬ 
mum  frequency  of  interest.  For  a  given  time  duration  the 
first  zero  crossing  occurs  at  the  lowest  frequency  for  a  square 
pulse.  For  that  type  of  pulse  the  first  zero  crossing  occurs 
at  a  frequency  equal  to  the  inverse  of  the  time  duration. 

A  good  rule  of  thumb,  then,  is  to  insure  that  the  duration 
of  the  impulse  is  less  than  2At,  where  t  is  the  sampling 
interval  in  the  analog-to-digital  conversion  process.  This 
would  put  the  first  zero  crossing  of  the  Fourier  transform 
of  a  square  pulse  at  the  Nyquist  folding  frequency,  and  the 
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first  zero  crossing  of  other  pulse  shapes  above  the  Nyquist 
folding  frequency. 

The  sample  length  is  equal  to  NAt  where  N  is  the 
number  of  digital  values  in  each  sample.  A  typical  value  of 
N  is  1024.  Thus,  the  duration  of  the  impulse  is  very  short 
relative  to  the  sample  length.  This  means  that  the  total 
energy  of  noise  represented  in  the  time-sample  can  be  on  the 
order  of  the  energy  of  the  impulse,  even  for  high  signal-to- 
noise  ratios.  The  noise  problem  is  further  aggravated  when 
employing  the  zoom  transform,  which  yields  increased  resolu¬ 
tion  in  a  given  frequency  band  by  effectively  increasing  the 
sample  length. 

With  other  techniques,  the  effects  of  noise  are 
reduced  by  averaging  the  power  spectrum  and  cross-spectrum 
functions  prior  to  the  computation  of  the  frequency  response 
function.  However,  only  a  few  averages  are  usually  used  in 
the  impulse  technique.  Otherwise,  the  time  advantage  of  the 
technique  is  lost.  Therefore,  special  time-sample  windows 
have  been  developed  for  the  impulse  technique. 

At  first  thought  it  might  seem  appropriate  to 
just  set  all  time-sample  values  beyond  the  impulse  to  zero, 
since  it  is  known  that  the  true  signal  value  after  the 
impulse  is  zero.  However,  this  would  be  equivalent  to  multi¬ 
plying  the  signal  by  a  narrow  rectangular  window.  In  applying 
any  type  of  window,  it  is  important  to  keep  in  mind  that 
multiplication  by  a  window  in  one  domain  is  equivalent  to 


convolution  of  the  Fourier  transforms  of  the  window  and  the 
data  in  the  other  domain,  resulting  in  distortion  of  the 
transformed  signal.  This  distortion  will  be  minimized  by 
minimizing  the  width  of  the  main  lobe  of  the  window  transform 
and  suppressing  its  side  lobes.  However,  there  is  a  funda¬ 
mental  conflict  between  these  requirements  and  the  reduction 
of  noise  in  the  time-sample  because  both  the  width  of  the 
main  lobe  and  the  amount  of  noise  reduction  are  inversely 
proportional  to  the  width  of  the  window  in  the  time  domain. 

To  further  complicate  the  situation,  suppression  of  the  side 
lobes  is  generally  achieved  at  the  expense  of  broadening  the 
main  lobe. 

A  good  compromise  has  been  arrived  at  in  practice 
in  the  form  of  a  window  with  unity  amplitude  for  the  duration 
of  the  impulse  and  a  cosine  taper,  with  a  duration  of  1/16 
of  the  sample  time,  from  unity  to  zero. 

Noise  problems  may  also  be  encountered  in  the 
response  signal,  particularly  when  dealing  with  heavily  damped 
systems  and  when  using  zoom  transform  analysis.  In  both 
cases  the  duration  of  the  response  signal  may  be  short  rela¬ 
tive  to  the  total  sample  time,  so  that  noise  may  comprise  a 
significant  portion  of  the  total  energy  in  the  time-sample 
even  with  relatively  high  signal-to-noise  ratios.  Another 
error  in  the  response  signal  that  is  encountered  when  testing 
lightly  damped  structures  occurs  when  the  response  signal  does 
not  significantly  decay  in  the  sample  window.  In  this  case 


the  resulting  time-sample  is  equivalent  to  multiplying  the 
true  response  signal  by  a  rectangular  window,  with  the  result 
that  the  frequency  resolution  may  not  be  sufficient  to  resolve 
individual  resonances. 

An  exponential  window  has  been  developed  to  re¬ 
duce  the  errors  that  occur  in  both  situations  described  above. 
The  window  decays  exponentially  from  1  to  a  value  of  0.05  in 
the  sample  time.  It  can  be  applied  directly  to  the  time- 
sample  of  the  response  signal  or  to  the  impulse  response 
function.  As  with  all  windows,  the  exponential  window  does 
change  the  resulting  measured  frequency  response  function; 
but  its  only  effect  is  to  increase  the  apparent  damping  in 
the  resonances.  It  does  not  change  the  resonance  frequencies 
and,  because  the  effect  of  the  exponential  window  is  the 
same  on  all  frequency  response  measurements,  it  will  not  alter 
the  measured  mode  shapes  if  applied  to  all  measured  frequency 
response  functions.  In  addition  to  reducing  noise  and  trun¬ 
cation  errors,  the  exponential  window  will  also  reduce  errors 
which  often  occur  when  testing  lightly  damped  systems  in 
which  the  damping  varies  with  the  measurement  position  on 
the  structure . 

Because  the  exponential  window  increases  the 
apparent  deamping  in  the  resonance  modes,  there  is  a  tendency 
of  the  window  to  couple  closely  spaced  resonance  modes.  Zoom 
transform  analysis  may  be  required  in  some  cases  to  allow 
sufficient  resolution  of  closely  spaced  modes  when  using  the 
exponential  window. 
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Zoom  transform  analysis  is  discussed  in  some  de¬ 
tail  along  with  several  examples  in  Reference  9.  It  is  a 
very  valuable  tool  in  impulse  testing,  as  it  is  in  other 
frequency  response  measurement  techniques.  The  effect  of  the 
zoom  transform  is  to  increase  the  resolution  of  the  analysis 
by  allowing  independent  selection  of  the  upper  and  lower 
frequency  limits  of  the  analysis  band.  With  the  zoom  trans¬ 
form,  for  example,  it  is  possible  to  perform  an  analysis  in 
the  frequency  range  from  900  to  1000  Hz  as  opposed  to  the 
corresponding  base-band  range  of  0  to  1000  Hz,  resulting  in 
a  10-to-l  increase  in  resolution,  for  a  given  sample  size  N, 
in  the  900  to  1000  Hz  band.  Because  of  greatly  increased 
resolution  possible  with  the  zoom  transform,  it  can  be  effec¬ 
tively  used  in  frequency  response  testing  to  separate  closely 
spaced  resonance  modes. 

There  are  two  important  effects  of  the  zoom 
transform  in  the  impulse  technique,  both  associated  with  the 
resulting  increase  in  sample  time.  The  first  effect  is  to 
make  possible  much  better  estimates  of  damping  in  lightly 
damped  systems.  This  is  due  to  the  reduction  of  truncation 
errors  in  the  sampled  response  signal.  The  second  effect, 
mentioned  previously,  is  aggravation  of  the  noise  problem 
in  both  the  input  and  response  signals.  The  second  effect 
makes  it  essential  that  force  and  response  windows  be  applied 
to  the  data  in  most  cases  when  using  the  zoom  transform  with 
the  impulse  technique. 


HP-54 51C  Fourier  Transfer  Function 


A  transfer  function,  as  determined  by  the  HP-5451C 
Fourier  Analyzer,  is  a  mathematical  description  of  a  system. 
It  can  be  defined  as: 


transfer  function 


Fourier  transform  of  output 
Fourier  transform  of  input 


or  equivalently. 


transfer  function 


average  cross  power  spectrum 

of  input  and  output _ 

average  power  spectrum 
of  input 


The  coherence  function  measures  the  degree  of  causality 
between  any  two  signals.  It  can,  therefore,  be  used  to  check 
the  validity  of  the  transfer  function.  When  a  transfer  func¬ 
tion  is  computed  it  may  not  be  obvious  that  there  are  ex¬ 
traneous  inputs,  or  that  the  system  is  nonlinear.  Both  of 
these  factors  would  introduce  error  in  the  computed  transfer 
function.  The  "Transfer  Function"  program  of  the  HP-5451C 
is  used  to  compute  the  transfer  and  coherence  functions.  A 
program  flow  chart  and  a  program  listing  are  provided  in 
Appendix  B. 


IV.  THEORY  OF  THE  TEST  PROCEDURE  TO  MEASURE  DAMPING 


The  design  concept  of  the  experimental  procedure  is  based 
on  the  Hewlett-Packard  modal  analysis  software  specifically 
designed  for  the  HP-5451C  Fourier  Analyzer.  The  theory 
of  the  complex  modes  for  damped  oscillatory  mechanical  sys¬ 
tems  is  described  below,  much  of  which  is  explained  in 
greater  detail  in  the  Modal  Analysis  Operating  and  Service 
Manual  {Option  402) . 

A.  MODAL  THEORY  OF  OPERATION  [REF.  6] 

Assume  that  the  motion  of  a  physical  system  can  be  des¬ 
cribed  by  a  set  of  n  simultaneous  second-order  linear  dif¬ 
ferential  equations  in  the  time  domain,  given  by, 

Mx  +  Cx  +  Kx  =  f  (9) 

where  the  dots  denote  differentiation  with  respect  to  time. 

f  =  f(t) 

is  the  applied  force  vector,  and 

x  =  x(t) 

is  the  resulting  displacement  vector,  while  M,  C,  and  K  are 
the  (n x n)  mass,  damping,  and  stiffness  matrices  respectively. 


Our  attention  will  be  limited  to  symmetric  matrices,  and  to 
real  element  values  in  M,  C,  and  K. 

Taking  the  Laplace  transform  of  the  system  equations  gives 

B (x)X(s)  =  Fla),  (10) 


where: 


B(S)  =  Ms2  +  Cs  +  K  (11) 

Here  s  is  the  Laplace  variable,  and  F(s)  is  the  applied  force 
vector  and  X(s)  is  the  resulting  displacement  vector  in  the 
Laplace  domain.  B(s)  is  called  the  system  matrix,  and  the 
transfer  matrix  H(s)  is  defined  as 

H  (s)  =  B  (s) _1  (12) 


which  implies  that 


H(s)F(s)  =  X(s)  (13) 

Each  element  of  the  transfer  matrix  is  a  transfer  function. 
The  elements  of  B  are  quadratic  functions  of  s,  and  since 
H  *  B-^,  it  follows  that  the  elements  of  H  are  rational 
fractions  in  s,  with  det(B)  as  the  denominator.  Thus,  H(s) 
can  always  be  represented  in  partial  fraction  form. 

If  it  is  assumed  that  the  poles  of  H,  i.e.,  the  roots  of 
det(B)  =  0,  are  of  unit  multiplicity,  then  H  can  be  expressed 
as 
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2n  ak 

H (s)  «  l  — g-  ,  (nxn)  (14) 

k=l  s  pk 

The  poles  occur  at  s  =  (zeros  of  det(B)),  and  each  pole 
has  n  (nxn)  residue  matrix  a^  associated  with  it.  For  an 
n*"*1  order  oscillatory  system,  there  will  always  be  2n  poles, 
but  they  will  appear  in  complex  conjugate  pairs.  Each  com¬ 
plex  pair  of  poles  causes  a  mode  of  vibration  in  the  struc¬ 
ture.  The  poles  are  complex  numbers  expressed  as 

pk  =  +  iuk  (15) 

where  ok  is  the  damping  coefficient  (a  negative  number  for 
stable  systems)  ,  and  is  the  natural  frequency  of  oscilla¬ 
tion.  The  resonant  frequency  is  given  by 


Q 


k 


(rad/sec) 


(16) 


and  the  damping  factor  is 


e 


(17) 


These  coordinates  are  shown  in  Figure  7.  When  5  =  1,  mode  (k) 
is  said  to  be  critically  damped.  It  is  also  possible  that 
5  >  1.  For  this  case  the  poles  of  mode  (k)  lie  along  the 
real  (or  damping)  axis  in  the  S-plane  and  it  is  said  to  be 
super-critically  damped. 
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Figure  7.  Poles  of  a  Mode  (k) 

Modal  vectors  (u^)  are  now  defined  as  solutions  to  the 
homogeneous  equation 


where : 


Bkuk 


=  0 


(10) 


Bk  =  B(pk) 

for  all  k  =  l,...,2n. 

In  addition,  pre-multiplying  B  times  the  equation  (14) 

and  letting  s  *  pk 

gives 


for  H(s),  multiplying  by  the  scalar  (s-pk) 
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(19) 


It  follows  by  comparison  of  (10)  with  (19)  that  every 
column  of  a^  must  contain  the  modal  vector  uk,  for  each 
k  =  1 ,  » » » , 2n . 

Similarly,  post  multiplying  B  times  the  equation  (14) 
for  H(s),  multiplying  by  the  scalar  (s-pk)  ,  and  letting 
s  =  pk  gives 

akBk  =0  (20) 

for  each  k  =  i,...,2n.  This  can  be  rewritten  as 

44  -  0  (21> 

where  t  denotes  the  transpose. 

But  Bk,  and  hence  ak,  is  assumed  to  be  symmetric  so 
equation  (21)  is  the  same  as  (19)  and  it  follows  that  each 
row  of  ak  must  also  contain  the  vector  uk  for  each  k  =  l,...,2n. 
In  order  to  satisfy  the  above  conditions  ak  must  have  the 

form 

aR  =  Akuku£  (n  x n)  (22) 

where  Ak  is  a  scalar. 

This  pervasiveness  of  the  modal  vectors  throughout  the 
transfer  function  matrix  is  evidence  of  the  so-called  global 
property  of  a  mode  of  vibration. 
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In  these  terms,  H  can  be  rewritten  as 


2n  A. 

H  =  V  — —  u,  uf 

kii  s_Pk  k  k 


(23) 


and  this  is  easily  written  in  matrix  form  as 


H  =  U  L  U‘ 


(n  x  n) 


(24) 


where  the  columns  of  U  comprise  the  u^  modal  vectors: 

U  *  [ux  u2  ...  u2n]  ,  (n  x  2n)  (25) 

and  L  is  a  diagonal  matrix  containing  all  s  dependence 


L  = 


Pre-mult iplying  H  by  Ut,  equation  (13)  ,  can  be  written  as 

(UfcUL) (UfcP)  =  (UfcX)  (27) 

so  that  transforms  the  spatial  vectors  F  and  X  to  vectors 

UfcF  and  UfcX  in  modal  coordinates.  Similarly  l^UL  is  the  modal 

representation  of  H.  Since  B(PJt)uk  *  0,  it  follows  that 
*  * 

B^Pk  uk  =  0'  so  t^xe  modal  vector  associated  with  the  conjugate 
pole  (p£)  is  u£  (the  conjugate  of  u^) . 

55 


s-p. 


2n 


s-p 


2n 


(2n  x  2n) 


(26) 


Thus,  the  above  U  matrix  always  contains  conjugate  pairs 
of  modal  vectors,  and  the  L  matrix  always  contains  elements 


corresponding  to  conjugate  pole  pairs  along  its  diagonal. 

If  is  defined  as  that  (n  *  n)  part  of  U  associated  with 
positive  poles,  then  will  correspond  to  the  negative  poles 
Similarly,  L  can  be  broken  into  two  parts,  comprising  the 
positive  poles,  and  Lj  comprising  the  negative  poles.  It 
can  then  be  represented 

H  =  U.L.U*  +  U*L2U*t  (28) 


or  in  partitioned  form  as 


H 


U 


*t 
1  ] 


(29) 


Each  of  these  sub-matrices  is  (n  x n)  and  only  and 
are  functions  of  s. 


A1 

s'pl 

• 

A1 

S-Pl 

• 

II 

r“t 

• 

• 

An 

S"Pn- 

n 

(M 

• 

• 

An 

s"pi 

(30) 


Then  H  can  be  written  as: 


H 


n 

l 

k=l 


u.uj;  J3L. 
k  k  s-pk 


*  *fc 

ukuk 


s-p. 


(31) 


Each  element  of  the  H  matrix  has  a  different  zero  in  the 
S-plane,  depending  upon  the  values  of  Ak  and  uk  at  each 
point,  but  the  poles  of  each  element  of  H  are  common,  and 
occur  at  s  =  pk  and  s  =  p*. 

1.  Identification  of  Modal  Parameters 

Because  of  the  form  of  the  ak  matrix,  only  one  row  or 
column  of  the  transfer  matrix  need  be  measured  and  analyzed, 
since  all  processes  of  measuring  the  transfer  matrix,  the 
unknown  parameters  in  equation  (14)  (i.e.,  the  complex  values 

of  pk  and  the  complex  values  of  the  elements  of  one  row  or 
column  or  the  residue  matrix  ak)  are  identified. 

Once  one  row  or  column  of  ak  has  been  identified,  it 
is  then  possible  to  construct  the  rest  of  the  rows  and  columns 

'f*  H 

in  a^.  For  example,  if  the  q  column  of  ak  is  given  by  a^, 
then 


a,  -  A,  u  ,  u, 
kq  k  qk  k 


(32) 


where  uqk  is  the  qtl1  component  of  the  modal  vector.  Since 
Akuqk  a  scalar»  clear  that  the  vector  of  modal 

residues  a^  is  proportional  to  the  modal  vector  u^,  for  each 
k  =  1, . . . ,2n. 


Since  A^  is  a  scaling  constant  it  can  be  assumed  that 
either  A^  *  1  or  the  square  root  of  u'u  =  1  without  loss  of 
generality  in  the  following  derivation. 
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Suppose  *  1,  then 


ak  =  ukuk 


and  for  the  q  column  or  row  of  a. 


akq  =  UqkUk 


Hence  the  q  element  of  av  is 

Kq 


akqq  ”  *Uqk^ 


Since  u  ^  is  a  scalar,  equation  (34)  can  be  rewritten 

akq 

uk  "  uT 
qk 

and  substituting  this  back  into  (33)  gives 


(Uqk> 


or  using  equation  (35) 


Hence  the  entire  matrix  a^  and  therefore,  the  entire 
transfer  function  matrix  H(s)  can  be  constructed  once  one  row 
or  column  of  residues  a^  has  been  identified  as  well  as  the 
pole  locations  pk,  for  each  k  =  l,...,2n.  The  residue  a. 
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is  at  the  driving  point  of  the  structure,  i.e.,  the  point 
where  the  structure  is  excited. 

2 .  Impulse  Response  of  Complex  Modes 

It  has  been  shown  that  a  mode  of  vibration  is  repre¬ 
sented  by  a  complex  conjugate  pair  of  poles  and  a  complex 
conjugate  pair  of  modal  vectors  in  the  transfer  matrix. 

Hence  for  a  single  mode  of  vibration  (k)  the  transfer  matrix 
is  written 


Hk(S) 


(n  x  n) 


(39) 


It  is  convenient  to  remove  a  factor  of  2i  from  the 
residue  matrix,  that  is  to  define  another  residue  matrix  rk 
such  that 


=  2. 


Then  the  transfer  matrix  is  written  as 


(40) 


Hk(s> 


K  K 

2i  (s-p,  )  "  2i  (s-p. ) 


(41) 


Each  component  of  this  matrix  exhibits  the  rectangular 
(or  coincident — quadrature)  from  shown  in  Figure  8  for  each 
valued  r^. 
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FREQUENCY 


REAL  PART 
H^UI 


HfcU) 


Figure  8.  Transfer  Function  for  a  Single  Mode  of 
Vibration 


Taking  the  inverse  Laplace  transform  of  (41)  gives 


x(t)  = 


rk  rx  X* 

e  -  »-r  e 


IT 


IT 


(42) 


-0.  t 

x(t)  =  |rk|e  sine (w^t  +  aR) 


(43) 


where  ak  »  the  angle  of  the  complex  residue  r^. 

Equation  (43)  is  the  impulse  response  of  a  complex 
mode  of  vibration.  The  mode  shape  is  defined  by  the  magnitude 

i  * 

I  9 

[  r.  ;  and  the  phase  angle  (av)  at  each  point  on  the  structure. 


The  mode  shape  matrix  is  printed  in  the  Modal  System  in 
terms  of  these  magnitude  and  phase  angles . 

3.  Modal  Mass,  Damping  and  Stiffness  and  Scaled  Mode 
Shapes 

Recall  that  the  assumed  symmetry  of  H(s)  along  with 
the  global  nature  of  mode  shapes  implies  that  the  residue 
matrix  for  each  mode  is  proportional  to  u^u^.  If  Qk  =  2iAk, 
then  equation  (22)  becomes 

rk  "  QkukUk  (nxn)  (44) 


Once  a  row  or  column  of  each  residue  matrix  r^  has 
been  identified,  the  mode  shape  vectors  can  be  calculated 
to  within  a  constant  of  proportionality.  This  vector  can 
be  scaled  large  or  small,  if  desired,  by  suitable  choice  of 
Qk*  Criteria  for  choosing  Qk  will  emerge  from  a  discussion 
of  modal  mass,  damping,  and  stiffness. 

Consider  a  single  degree  of  freedom  system  represented 
by 


mx  +  cx  +  kx  *  f 


(45) 


where  m,  c,  and  k  are  scalars.  The  Laplace  Transform  gives 

(ms^  +  cs  +  k)X(s)  =  F(s)  (46) 


so  that 


H  (s) 


X(s) 

FlxT 


(47) 
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If  the  damping  factor  is  less  than  one,  the  polynomial 


in  the  denominator  can  be  factored. 


H(s) 


1 

_ m 

(s-p) (s-p*) 


where 


_C  _  v//  C  .  2  _  k 
2m  '2m'  m 


(48) 


(49) 


P 


-c 

2m 


+  i 


m 


(  c  )2 
(5m} 


(50) 


p  =  -o  +  ini,  o  >  0,  u>  >  0  (51) 


and 


* 

p  =  -o  -  iu) 


(52) 


To  express  Equation  (48)  in  the  partial  fraction  form 
of  Equation  (39),  solve  for  a  and  a*. 


H  (s) 


1 

_ m _ 

(s-p) (s-p*) 


a 

(s-p) 


+ 


* 

a 

(S-p*) 

(53) 


1 

m 

(s-p*) 


3— P 


[a  + 


a* (s-p) , 
*  J 

(s-p  ) 


s=p 


(54) 
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1.: 


1 


a 


(55) 


1 


m 

-2iw 


★ 


a 


(56) 


H  (s) 


1 


mui 

2i(s-p) 


+ 


JL_ 

mm 

2i (s-p*) 


(57) 


The  residue  r  can  be  found  from  Equations  (40),  (41), 


and  (57). 


r  = 


mw 


(58) 


Now  note  that  Equations  (50)  ,  (51) ,  and  (58)  give  m, 
c,  and  k  in  terms  of  r,  u,  and  o. 


m  = 


rui 


(59) 


c  =  2  am 


(60) 


2  2 

k  -  (ai  +  a  )m 


(61) 


For  multiple  degree  of  freedom  systems,  the  analogous 
definitions  are  made  for  each  mode  k. 


m 


k 


1 

Vk 


(62) 


2°kmk 


(63) 
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(64) 


kk  =  +  °k)mk 

Equations  (44)  ,  (62)  ,  (63)  ,  and  (64)  are  used  to 
calculate  modal  mass,  damping,  and  stiffness  and  scaled  mode 
shapes  in  the  Print  Step.  Qk  is  arbitrary  so  it  can  be 
chosen  to  give  unit  mass,  unit  u*"u,  etc. 

The  modal  mass,  stiffness,  and  damping  coefficients 
can  be  interpreted  as  being  the  masses,  dampers,  and  springs 
of  decoupled,  single  degree  of  freedom  systems.  These  systems 
are  equivalent  to  the  original  system  under  a  change  of 
coordinates . 

4 .  Measurement  Implications  of  Modal  Theory 

A  fundamental  assumption  of  modal  testing  is  that  a 
mode  of  vibration  can  be  excited  from  anywhere  on  an  elastic 
structure,  except  of  course  along  its  node  lines  (zero  points) 
where  it  can't  be  excited  at  all.  This  is  another  way  of 
stating  the  result  derived  earlier  (i.e.,  that  the  same  modal 
vector — scaled  by  a  different  component  of  itself--is  con¬ 
tained  in  every  row  and  column  of  the  transfer  matrix) .  In 
addition,  modal  frequency  and  damping  are  constants  which  can 
be  identified  in  any  element  of  the  transfer  matrix,  i.e., 
any  transfer  function  taken  from  the  structure . 

It  is  important  to  recognize  that  this  global  mode 
shape  concept  exists  within  certain  spatial  boundaries,  beyond 
which  vibrations  will  not  readily  propagate.  If  two  linear 
systems  are  completely  isolated,  then  a  single  composite  mode 
including  both  systems  is  not  meaningful. 
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Conversely,  it  is  important  to  include  enough  spatial 
points  in  the  measured  data  set  to  describe  all  of  the  vibra¬ 
tion  modes  of  interest.  If  some  region  of  a  bounded  system 
is  not  monitoried  or  excited,  or  if  points  are  not  chosen 
sufficiently  close  together,  then  some  modes  cannot  be  ade¬ 
quately  represented. 

When  a  system  is  represented  by  the  partial  fraction 
form  of  its  transfer  matrix,  a  closed  form  solution  for  the 
displacement  at  any  point,  for  any  combination  of  modes,  is 
readily  obtainable  using  simple  matrix-vector  multiplication. 
This  is  particularly  helpful  when  the  response  is  only  to  a 
few  modes  of  interest. 

The  primary  purpose  for  using  modal  coordinates  to 
describe  the  dynamics  of  a  linear  vibration  structure  is  that 
it  drastically  reduces  the  amount  of  time  and  effort  neces¬ 
sary  to  measure  the  dynamics  in  a  laboratory.  The  essence  of 
the  modal  concept  is  that  once  one  row  or  column  of  the 
transfer  function  matrix  has  been  determined,  the  entire  matrix 
and  hence  the  entire  dynamics  of  the  structure  can  be  specified. 

In  order  to  obtain  valid  modal  results  with  the  Modal 
Analysis  system,  several  important  assumptions  of  the  modal 
theory  described  here  must  be  satisfied. 

1.  The  structure  must  exhibit  the  behavior  of  a  linear 
system.  If  you  are  not  able  to  successfully  curve 
fit  equation  (14)  to  measured  transfer  function 
data,  then  the  behavior  may  not  be  linear. 
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The  structure  must  exhibit  the  reciprocity  of  symmetry 
property.  This  can  be  verified  by  comparing  the 
transfer  function  obtained  from  an  excitation  measure¬ 
ment  at  point  A  and  response  measurement  at  point  B 
with  the  transfer  function  obtained  from  an  excitation 
at  point  B  and  response  at  point  A.  This  should  be 
the  same. 

The  structure  must  have  distinct  pole  pairs.  Since 
the  Modal  Analysis  system  uses  analytical  expression 
(6)  to  perform  the  curve  fitting,  it  cannot  fit  re¬ 
peated  poles  (repeated  roots  of  det(B)  =  0  because 
this  expression  was  derived  under  the  assumption  that 
the  poles  are  unique.  This  condition  can  be  difficult 
to  detect  in  one  set  of  transfer  function  measurements. 
To  insure  that  the  poles  are  unique  you  should  measure 
a  different  row  or  column  of  the  transfer  function 
matrix,  and  compare  the  mode  shapes  from  the  two  sets 
of  measurements.  If  they  are  not  identical  there  is 
a  strong  possibility  of  repeated  roots. 

There  may  be  more  than  one  modal  vector  corresponding 
to  a  given  pole,  due  to  certain  kinds  of  boundary 
conditions  on  the  structure,  i.e.,  symmetric  conditions 
along  several  axes.  To  insure  that  this  is  not 
happening,  you  should  compare  modal  vectors  obtained 
from  two  or  more  rows  or  columns  of  the  transfer 


matrix . 


5.  The  transfer  function  H(s)  is  defined  in  terms  of  an 


input  force  and  a  response  displacement.  Response 
velocity  or  acceleration  can  be  converted  to  displace¬ 
ment.  But  non- force  Modal  Analysis  system  will  calcu¬ 
late,  print,  and  display  modal  data  (i.e.,  poles  and 
residues)  for  such  measurements,  however.  The  excep¬ 
tions  are  the  four  mode  shape  scaling  procedures  that 
produce  modal  mass  damping  and  stiffness.  The  system 
logica  will  reject  an  attempt  to  calculate  these  values, 
unless  motion  in  response  to  force  was  measured. 
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4 


r.  PROCEDURE  OF  TEST  CHAMBER  CHARACTERIZATION 

A.  PROCEDURE  OF  DATA  COLLECTION 

After  construction  of  the  test  chamber,  preservation  and 
painting,  installation  of  cooling/heating  coils,  and  instal¬ 
lation  of  all  insulation,  the  following  is  the  procedure  used 
for  final  characterization  of  the  test  chamber  (with  front 
panel  removed) . 

1.  Draw  a  2.54  inch  *  2.54  inch  grid  on  entire  inside 
surface  of  the  test  chamber.  This  grid  is  used  to  record 
location  of  impact  or  excitation,  as  well  as  location  of 
response  transducer  (Figure  9) . 

2.  Re-install  exciter  and  support  framework;  re-install 
specimen  fixture  with  specimen  in  place;  and  re-install 
test  chamber  top.  Drill  1/4  inch  hole  through  specimen  for 
attachment  of  exciter  driver  rod.  Install  exciter  driver 
rod  with  Force  Transducer  in  place  (Figure  10) .  NOTE: 
Acoustic  exciter  and  force  transducer  will  not  be  used  in 
this  work.  For  impulse  hammer  testing  disconnect  exciter 
from  specimen. 

3.  Using  accelerometer  mounting  wax,  attach  pick-up  trans¬ 
ducer  to  desired  location.  Connect  signal  conditioners  and 
wiring  as  instructed  in  the  operating  manuals.  System  flow 
chart  is  presented  in  Figure  11.  Impulse  hammer,  signal 
conditioners  and  typical  attachment  of  pickup  transducer  are 


fntce 


transducer 


Figure 


10. 
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IMPULSE  RESPONSE 


Figure  11.  Impulse  hammer  technique  flow  chart 
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presented  in  Figures  12  and  13.  All  calibration  data 
of  sensors  and  mounting  and  instructions  are  included 
in  Appendix  C. 

4.  Turn  on  the  HP-5451C  Fourier  Analyzer  (HP-2648A  Graphics 
Terminal  will  automatically  be  energized) .  Insert 
"Modal"  disk  into  disk  drive.  Energize  disk  drive. 

When  "Disk  Ready"  light  comes  on  "Boot"  the  system  as 
outlined  in  volume  1  of  System  Operating  Manual  (Figures 
14  and  15) .  Activate  user  assignable  keys  fl  and  f2 

on  the  graphics  terminal  to  "Graphics  Mode"  and  "Alpha 
Mode"  respectively.  Procedures  for  key  assignments  are 
found  in  volume  5  of  the  Computer  Operating  Manual. 

5.  From  this  point  on,  the  procedure  for  taking  impulse 
data  is  interactive.  Locally  generated  computer  pro¬ 
grams,  which  are  designed  to  allow  the  data  generated 
by  impulse  hammer  technique  to  be  stored  and  later 
processed  by  the  off-line  Modal  software,  are  presented 
in  Appendix  D.  What  follows  is  a  typical  step-by-step 
impulse  data  taking  session  assuming  pick-up  transducer 
at  location  L051  and  hammer  impulse  location  at  L066 
The  following  abbreviations  will  be  used: 

C  =  computer  keyboard  entry 

D  =  graphics  terminal  display 

G  =  graphics  terminal  entry 

R  =  user  response 


72 
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HP-54420A  LOW-PASS 
FILTERS  AND  ANALOG- 
DIGITAL  CONVERTER 


•HP-1000  COMPUTER 


IP-5477A  MODAL  SECTION 


.HP5460A  DISPLAY  SECTION 


.COMPUTER  KEYBOARD 


< 


HP-54 66B  DIGITAL- ANALOG 
CONVERTER 


.HP-7960  DISK  DRIVE 


5451C  Fourier  Analyzer 
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Step  Number 


1. 


2. 


3. 


4. 

5. 


6. 


7. 


8. 


9. 


10. 


11. 


12. 


13. 


14. 


15. 


Normal  "Boot”  of  modal  with  fl  &  f2  programmed. 

D:  SET  UP  DATA  NUMBER?  (1  to  5) . 

R;C:  2  [ENTER]  (want  to  preserve  Demonstration  data 

in  set  up  #1,  set  ups  #2  through  #5  are  available) 

D:  LIST (1 ) ,  EDIT(2),  FINISH(3)? 

R;C:  2  [ENTER]  (must  enter  particular  transducer 
specifications  and  general  information  for 
a  particular  test) . 

D:  LABEL? 

R;G:  TANK  CHARACTERIZATION  [RETURN]  (title  of 

particular  test)  . 

D:  LINE  NOS.?  (FIRST , LAST) . 

R;G:  1,6  [RETURN]  (you  may  modify  one  or  any 

sequential  series  of  lines  of  the  set  up  data) . 

D:  1.  NO.  OF  MEASUREMENTS? 

R;G:  16  [RETURN]  (A  measurement  is  one  record  of 

storage  on  the  disk.  One  record  is  1024  words 
in  length,  and  a  test  is  n  measurements  stored 
in  n  contiguous  records  on  the  disk) . 

D:  STARTING  LOCATION  OF  MEAS .  STORAGE  (1  to  485). 

R;G:  5  [RETURN]  (storage  location  1  through  255  are 

on  the  removable  disk,  and  storage  locations 
266  to  485  are  on  the  fixed  disk) . 

D:  DATA  BLOCK  SIZE?  (64  to  2048) . 

R;G:  1024  [RETURN]  (block  size  of  1024  will  give 

maximum  zoom  capability) . 
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16.  D:  4.  POLAR  TRANS.  FUNCT.  DISPLAY?  LINEAR (1)  or 

LOG (2)  MAGNITUDE. 

17.  R;G:  2  [RETURN]. 

18.  D:  5.  TRANSDUCER  UNITS?  ENGLISH (1) ,  METRIC(2). 

19.  R;G:  1  [RETURN]. 

20.  D:  TRANSDUCER  SCALE  FACTORS  (VALUE, UNITS) 

1  *  MV/G,  2  =  VM/IPS,  3  =  MV/IN,  4  =  MV/LB 
A.  INPUT? 

21.  R;G:  10.6,4  [RETURN]  (values  are  obtained  from 

calibration  data  of  sensors) . 

22.  D:  B.  RESPONSE? 

23.  R;G:  10.12,4  [RETURN]. 

24.  D:  G.  AMPLIFIER  GAINS? 

A.  INPUT? 

25.  R; G :  1  [RETURN]  (enter  transducer  amplifier  switch 

setting) . 

26.  D:  B.  RESPONSE? 

27.  R;G:  1  [RETURN] 

28.  D:  LIST ( 1) ,  EDIT (2),  FINISH(3)? 

NOTE:  At  this  point  (1)  may  be  entered  and  all  data  may 
be  checked.  If  any  corrections  are  required,  enter  (2) 
and  use  procedures  outlined  in  steps  5  through  27. 

29.  R;G:  3  [RETURN]. 

NOTE:  The  "set  up"  portion  of  the  test  is  now  complete. 
The  next  step  is  the  "measure"  portion. 

30.  D:  STEP  NUMBER?  (1  to  5). 


31.  R;C:  2  [ENTER]  (the  computer  is  now  in  the  "measure" 

mode) . 

32.  D:  TRANS ( 1)  ,  RANDOM ( 2 )  ,  DISK(3)/  USER(4), 

FOURIER/ZOOM/GRAPHICS (5) ? 

33.  R;C:  5  [ENTER]  (the  HP  5451-C  is  now  ready  to  set 

up  for  an  impulse  test) . 

34.  R;C:  SELECT  BLOCK  SIZE  (1024  MAX  FOR  ZOOM) 

SET  ADC  FREQUENCY  RANGE  AS  DESIRED 

SET  TREGGER  TO  INTERNAL 

35.  R;C:  [JUMP]  0  [SPACE]  1  [ENTER]  (this  command  calls 
the  user  programs  1,50,51,52,  and  59  for  use). 
TRANSFER  FUNCTION  MEASUREMENTS 

ARE  HP  FILTERS  INSTALLED? 

36.  D:  (0  -  NO  1  -  YES) 

37.  R; C:  1  [ENTER] 

38.  D:  CUT-OFF  FREQUENCY  (CHANNEL  NO)? 

NOTE:  To  eliminate  unwanted  or  meaningless  data  beyond 

range  of  impulse  hammer,  enter  channel  number  above  which 
the  computer  will  not  sampel  data. 

EXAMPLE:  0-500  lb.  hammer  useful  range  is  from 

0-6000  Hz,  with  block  size  1024  selected 
channels  above  307  are  meaningless. 

39.  R; C:  307  [ENTER] 

40.  D:  IMPACT  (1)  OR  RANDOM  (2)  EXCITATION? 

41.  R; C:  1  [ENTER] 

42.  D:  NO.  OF  AVERAGES? 


43.  R; C:  4  [ENTER] 

44.  D:  BASEBAND  OR  ZOOM  MEASUREMENT? 

(0  =  BASEBAND  1  =  ZOOM) 

45.  R;C:  0  [ENTER]  (first  measurement  must  be  baseband). 

46.  D:  BASEBAND  IMPACT  MEASUREMENT 
PRESS  CONTINUE  FOR  MEASUREMENT 

NOTE:  At  this  point  trigger  and  analog-Digital  Converter 
levels  need  to  be  set  to  do  this.  The  following  keg 
strodes  should  be  done. 

47.  R;C:  SET  REPEAT/SINGLE  SWITCH  TO  REPEAT 
PRESS  ANALOG  IN 

SET  ADC  RANGE  SWITCHES  AND  TRIGGER  SLOPE  AND  LEVEL 
SET  ADC  TO  SINGLE  WHEN  READY  AND  PRESS  CONTINUE 

48.  R; C:  [CONTINUE]. 

NOTE:  Impact  4  times.  The  computer  will  beep  after  each 

data  entry.  After  the  computer  processes  the  data  the 
ready  light  will  come  on  and  the  data  in  Block  0  will 
be  displayed  on  the  small  screen  of  the  computer. 
Processed  data  are  stored  in  the  following  locations  and 
can  be  displayed  by  the  key  strokes: 

(DISPLAY] 0  [ENTER]  Log  transfer  function 
[DISPLAY] 1  [ENTER]  Coherence 
[DISPLAY] 2  [ENTER]  Input  power  spectrum 
[DISPLAY] 3  [ENTER]  Output  power  spectrum 
[DISPLAY] 4  [ENTER]  Cross  power  spectrum 
[DISPLAY] 5  [ENTER]  Raw  data 


To  obtain  hard  copy  of  any  of  the  above  outputs  execute  the 
following  key  commands  [USERl  5821  [SPACE]  35  [ENTER] .  Then 
enter  [USER]  [PLOT]  [ENTER] .  To  return  plotting  to  graphics 
execute  [USER]  5821  [SPACE]  6  [ENTER] ,  and  when  the  Graphics 
Terminal  is  in  the  Graphics  mode  (depress  fl)  and  [USER] 

[PLOT]  [ENTER]  is  executed  all  plotting  to  be  sent  to  the 
graphics  terminal.  To  continue,  f2  must  be  depressed.  After 
all  preliminary  graphing  is  completed  on  baseband  data  press 
[continue] .  Typical  plots  of  baseband  data  are  presented  in 
Figures  16  through  20. 

49.  D:  PRESS  CONTINUE  WHEN  READY 

50.  R;C:  [CONTINUE] 

51.  D:  SAVE  DATA? 

(0  =  NO  1  =  YES) 

52.  R; C:  0 [ENTER] 

53.  D;  MAKE  ANOTHER  MEASUREMENT? 

()  =  NO  1  =  YES)  . 

54.  R; C:  1 [ENTER] 

55.  D:  NO.  OF  AVERAGES? 

56.  R; C :  R [ENTER] 

57.  D:  BASEBAND  OR  ZOOM  MEASUREMENT? 

()  =  BASEBAND  1  =  ZOOM) 

58.  R; C:  1 [ENTER] 

59.  D:  ZOOM  IMPACT  MEASUREMENT 

MOVE  CURSOR  TO  START  FREQUENCY 
PRESS  "VALUE"  (SWITCH  REGISTER  11) 
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power  spectrum  of  baseband  data 


V  SQR 


Figure  20.  Cross  power  spectrum  of  baseband  data 


NOTE:  The  desired  start  and  end  frequency  for  the  zoom 
measurement  is  not  entered  with  use  of  cursor. 


60.  D:  "TYPICAL"  ZOOM  IMPACT  MEASUREMENT 

MOVE  CURSOR  TO  START  FREQUENCY 
PRESS  "VALUE"  (SWITCH  REGISTER  11) 
CHANNEL  =  49.0000 

FREQUENCY  =  957.0999  HZ 
AMPLITUDE  =  -700.0026  E  -9.0000 
MOVE  CURSOR  TO  END  FREQUENCY 
PRESS  "VALUE" 

61.  D:  CHANNEL  =  134.0000 

FREQUENCY  =  2618.0000  HZ 

AMPLITUDE  =  -381.7491  E  -9.0000 

ANALYZE  OLD  OR  NEW  DATA? 

1  =  OLD  (FROM  THROUGHPUT  FILE) 

2  =  NEW  (WILL  WRITE  TO  THROUGHPUT  FILE) 

PRESS  CONTINUE  FOR  MEASUREMENT 

62.  R; C:  2  [ENTER] 

NOTE:  Whenever  zoom  is  being  done  new  data  is  required. 

63.  R;C:  [CONTINUE] 

NOTE:  For  zoom  measurement  the  impact  should  be  a  series 
of  rapid  impacts  with  a  varying  interval  between  impacts. 
The  trigger  light  will  remain  lit  thorughout  the  data 
taking  session. 

64.  D:  ANALYZING  THROUGHPUT  DATA  "TYPICAL" 

CNTR  FREQ  1479  HZ/DIV  125.0 
DF:  000.2441406  BLOCKS  LEFT  25 


ZOOM  POWER  4 


NOTE:  After  data  has  been  analyzed  and  desired  plots 
made  as  described  in  step  43,  ZOOM  data  should  be  stored 
for  later  use  by  the  off-line  modal  software.  A  rec- 
commended  data  sheet  is  shown  in  Figure  21.  It  is 
important  to  log  number  of  measurements  and  start  loca¬ 
tion  of  measurement  storage. 

65.  D:  PRESS  CONTINUE  WHEN  READY 

NOTE:  Generally  it  is  desired  to  make  multiple  tests  at 
different  locations  at  the  same  frequency  range  so  that  a 
modal  analysis  can  be  done  at  a  later  time. 

66.  R;C:  [CONTINUE] 

67.  D:  SAVE  DATA? 

(0  =  No  1«  YES) 

NOTE:  The  data  taking  session  can  be  repeated  by  entering 
"0" . 

68.  R;C:  1 [ENTER] 

69.  D:  DATA  STORAGE  NO.? 

70.  R;C:  XXX [ENTER] 

NOTE:  Data  storage  registers  0  through  500  are  available 

for  permanent  storage  of  the  data. 

71.  D:  MAKE  ANOTHER  MEASUREMENT? 

(0  =  NO  1  =  YES) 

NOTE:  From  this  point  the  steps  57  through  71  will  be 
repeated  for  as  many  zoom  measurements  as  desired.  Each 
zoom  measurement  will  be  stored  in  a  next  higher  record 
number  with  the  begin  record  number  known. 
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After  all  data  has  been  taken  and  recorded,  and  the 
number  of  measurements  and  start  record  measurement  has 
been  correctly  entered  in  the  set-up  table,  the  system  is 
ready  to  proceed. 

72.  R;C:  [SELECT  STEP]  Modal  Analysis  Controller 

73.  D:  STEP  NUMBER? 

74.  R; C :  [X]  [ENTER]  where  X  =  set  up 

data  number  i.e., 

1  for  example 

75.  D:  LIST (1) ,  EDIT(2),  FINISH(3)? 

76.  R;C:  3 [ENTER]  assuming  all  data  correct. 

From  this  point  on,  the  procedure  for  the  mode  identifi¬ 
cation  and  damping  measurement  is  covered  in  Reference  6 . 

For  the  tank  characterization  the  modes  identified  and  damping 
measured  at  the  particular  modes  are  presented  in  Table  II. 


B.  TEST  CHAMBER  MAJOR  MODES  AND  DAMPING 

To  identify  the  major  modes  and  measure  the  damping  of 
the  test  chamber  the  polar  plot  of  the  baseband  band  was 
divided  into  10  sections  (Figure  22) ,  and  each  section  was 
zoomed  using  the  same  pickup  location  and  using  the  same  im¬ 
pact  location  on  left,  right,  deck  and  back  panel. 
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of  the  test  chamber 


TABLE  I 


TEST  CHAMBER  FREQUENCY 

SECTIONS 

SECTION 

START  FREQ. 

STOP  FREQ. 

CNTR.  FREQ. 

DELTA  FREQ 

NO. 

(Hz) 

(Hz) 

(Hz) 

(Hz) 

1 

423 

1,138 

781 

1.395 

2 

976 

1,471 

1,221 

0.977 

3 

1,583 

2,247 

1,780 

1.786 

4 

2,247 

2,715 

2,481 

1.221 

5 

2,735 

3,321 

3,028 

1.627 

6 

3,340 

3,711 

3,525 

0.977 

7 

3,711 

4,493 

4,102 

1.953 

8 

4,473 

4,805 

4,639 

0.814 

9 

4,805 

5,411 

5,108 

1.630 

10 

5,411 

5,723 

5,567 

0.814 

A  modal  analysis 

of  each  section 

was  completed. 

Major 

modes  and  damping  for 

each  major  mode 

was  analyzed. 

The 

results  are  presented  in  Table  II.  Test  data  are  included 
in  Appendix  E. 
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TABLE  II 


TEST  CHAMBER  MODES  AND  DAMPING  FACTORS 


SECTION 

NO. 

NO.  OF 

MODES 

AVERAGE  DAMPING 
FACTORS  (%) 

STANDARD  DEVIATION 
OF  DAMPING 

1 

3 

0.5895 

0.3059 

2 

2 

0.3642 

0.0035 

3 

10  ' 

0.8918 

0.2458 

4 

12 

0.7375 

0.5084 

5 

10 

0.3601 

0.1934 

6 

6 

0.2763 

0.1065 

7 

8 

0.3115 

0.2176 

8 

3 

0.2219 

0.0627 

9 

10 

0.1787 

0.0873 

10 

2 

0.1660 

0.0525 

The 

overall  average 

damping  factor  of 

the 

test  chamber  is 

0.4667% 

with  an  overall 

standard  deviation 

of 

0.3689%. 
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VI.  PROCEDURE  FOR  DAMPING  MEASUREMENT  OF  SPECIMEN 

A.  SPECIMEN  SECURED  IN  FIXTURE  INSIDE  TEST  CHAMBER 

The  test  specimen  (cast  nickel-aluminum  bronze,  code 
FTC)  was  secured  in  the  test  chamber,  impact  locations  A,  B, 

C,  are  shown  in  Figure  23.  To  identify  the  major  modes  and 
measure  the  damping  factor  of  the  specimen  the  0  to  500  lb. 
impulse  hammer  was  used  for  the  frequency  range  up  to  approxi¬ 
mately  6000  Hz,  and  the  0  to  50  lb.  impulse  hammer  was  used 
for  the  frequency  range  from  approximately  5000  Hz  to  approxi¬ 
mately  12,000  Hz.  Baseband  polar  plots  for  each  impulse 
hammer  at  each  location  are  shown  in  Figures  24  through  29. 

It  should  be  noted  that  there  is  very  good  correlation  between 
the  0-500  lb.  impulse  hammer  and  the  0-50  lb.  impulse  hammer. 
Additionally,  there  is  very  good  correlation  between  various 
impact  locations. 

Utilizing  the  baseband  data  (0-500  lb.  impulse  hammer), 
obtained  for  impact  location  B,  five  frequency  ranges  were 
investigated.  Additionally,  using  the  baseband  data  obtained 
from  the  0-50  lb.  impulse  hammer  at  location  B,  the  frequency 
ranges  of  three  additional  ranges  were  investigated.  Specimen 
Frequency  sections  are  listed  in  Table  III.  Specimen  modes 
and  damping  factors  are  presented  in  Table  IV.  Measurement 
data  and  rough  log  sheets  are  included  in  Appendix  F. 
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Cast  nickel-aluminum  bronze,  code  FTC  specimen 
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Figure  24.  Polar  presentation  of  transfer  function  for  baseband 
(0-500  lb.  impulse  hammer).  Impact  location  A 


Figure  25.  Polar  presentation  of  transfer  function  for  baseband 
(0-50  lb.  impulse  hammer) .  Impact  location  A 


Figure  26.  Polar  presentation  of  transfer  function  for  baseband 
measurement  (0-500  lb.  impulse  hammer).  Impact 
location  B.  Zoom  ranges  identified 


Figure  27.  Polar  presentation  of  transfer  function  for  baseband  measurement 

(0-50  lb.  impulse  hammer) .  Impact  location  B.  Zoom  ranges  identified 


Figure  29.  Polar  presentation  of  transfer  function  for  baseband  measurement 
(0-50  lb.  impulse  hammer) .  Impact  location  C 


TABLE  III 


SPECIMEN  FREQUENCY  SECTIONS  (SPECIMEN  FIXED) 


SECTION 

NO. 

START  FREQ. 
(Hz) 

STOP  FREQ. 
(Hz) 

CNTR.  FREQ 
(Hz) 

DELTA  FREQ 
(Hz) 

1 

410 

937 

673 

1.395 

2 

840 

2,598 

1,758 

4.882 

3 

2,598 

3,340 

2,969 

1.953 

4 

4,043 

4,512 

4,277 

1.221 

5 

4,747 

5,508 

5,127 

1.953 

6 

5,713 

7,276 

6,494 

4.069 

7 

7,276 

9,620 

8,448 

6.104 

8 

9,620 

11,380 

10,500 

4.883 

Modal  analysis  of 

the  results  are  shown 

the  specimen  sections  was 

in  the  following  table 

completed  and 

TABLE  IV 


SPECIMEN  MODES  AND  DAMPING  FACTORS  (SPECIMEN  FIXED) 

SECTION  NO.  OF  AVERAGE  DAMPING  STANDARD  DEVIATION 
NO.  MODES  FACTORS  ( % )  OF  DAMPING  FACTOR 

1  1  0.1755  - 

2  2  0.2173  0.1088 

3  3  0.0559  0.0063 

4  2  0.0883  0.0025 

5  2  0.0958  0.0576 

6  6  0.1310  0.0542 

7  5  0.0917  0.0541 

8  5  0.1000  0.0396 

The  overall  average  damping  factor  for  the  specimen  is 
about  0.1112%  with  an  overall  standard  deviation  of  about 


0.0601% 


B.  SPECIMEN  REMOVED  FROM  FIXTURE  AND  TEST  CHAMBER 

In  order  to  ascertain  the  effects  of  the  specimen  fix¬ 
ture  and  test  chamber,  one  additional  test  was  conducted  on 
the  same  specimen.  During  this  test  the  specimen  was  removed 
from  the  fixture  and  test  chamber.  The  specimen  was  laid 
flat  on  3/4  inch  foam  rubber.  No  other  support  was  provided. 

Response  and  impulse  hammer  locations  are  identical  to 
the  previous  test.  Frequency  ranges  for  use  of  the  large 
(0-500  lb.)  impulse  hammer  and  small  (0-50  lb.)  impulse  hammer 
remained  unchanged. 

The  baseband  plots  for  large  and  small  impulse  hammer 
for  locations  A,  B,  and  C  are  presented  in  Figures  30  through 
35.  For  the  large  impulse  hammer  four  ranges  of  zoom  were 
completed,  and  for  the  small  impulse  hammer  five  additional 
sections  were  zoomed.  Specimen  frequency  sections  are  tested 
in  Table  V. 
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Figure  31.  Polar  presentation  of  transfer  function  for  baseband 

measurement  (0-50  lb.  impulse  hammer) .  Impact  location 
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Figure  32.  Polar  presentation  of  transfer  function  for  baseband  measurement 

(0-500  lb.  impulse  hammer) .  Impact  location  B.  Zoom  ranges  identified 


.  -  iij.  iiij 


Figure  33.  Polar  presentation  of  transfer  function  for  baseband  measurement 

(0-50  lb.  impulse  hammer).  Impact  location  B.  Zoom  ranges  identified 


Figure  34.  Polar  presentation  of  transfer  function  for  baseband  measurement 
(0-500  lb.  impulse  hammer).  Impact  location  C 


Polar  presentation  of  transfer  function  for  baseband  measure' 
ment  (0-50  lb.  impulse  hammer) .  Impact  location  C 


TABLE  V 


SPECIMEN  FREQUENCY  SECTIONS  (SPECIMEN  FREE) 


SECTION 

NO. 

START  FREQ. 
(Hz) 

STOP  FREQ. 
(Hz) 

CNTR  FREQ. 
(Hz) 

DELTA  FREQ. 
(Hz) 

1 

1094 

1954 

1524 

2.441 

2 

1915 

3555 

2735 

4.883 

3 

3575 

4610 

4092 

3.255 

4 

4610 

5293 

4951 

1.953 

5 

4688 

6299 

5493 

4.069 

6 

6250 

7227 

6738 

2.441 

7 

7276 

8594 

7935 

3.441 

8 

8545 

9473 

9009 

2.441 

9 

9424 

11,190 

10,307 

4.883 

Modal  analysis  was  compled  on  the  free  specimen,  and  the 
results  are  tabulated  below  in  Table  VI. 


( 
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TABLE  VI 


SPECIMEN  MODES  AND  DAMPING  FACTORS  (SPECIMEN  FREE) 


SECTION  NO.  OF  AVERAGE  DAMPING  STANDARD  DEVIATION 
NO.  MODES  FACTORS  (%)  OF  DAMPING  FACTOR 


0.0885 

0.1264 

0.0339 

0.0265 

0.0295 

0.0250 

0.0331 

0.0209 

0.0301 


0.0064 

0.0774 

N/A 

0.0017 

0.0020 

0.0080 

0.0195 

0.0123 

0.0069 


Again,  all  raw  data  sheets  are  included  in  Appendix  F. 
The  overall  average  damping  factor  for  the  free  specimen  is 
about  0.0500%  with  an  overall  standard  deviation  of  about 
0.0504%. 


VII.  RESULTS  AND  CONCLUSIONS 


Mode  identification  was  accomplished  utilizing  the  zoom 
feature  of  the  HP-5451C  Fourier  Analyzer.  After  a  baseband 
test  was  conducted  and  results  stored  (Figure  22) ,  any  section 
could  be  chosen  to  zoom.  For  example,  zooming  section  1  of 
the  test  chamber  results  in  a  plot  as  shown  in  Figure  36. 

The  initial  display  will  be  in  polar  form  and  must  be  con¬ 
verted  to  rectangular  (one  button  function  on  the  HP-5451C) 
as  shown  in  Figure  37.  As  shown  in  Figure  38,  the  coherence 
for  this  test  was  very  good,  and  the  data  may  be  assumed  to 
be  accurate.  For  mode  identification  it  is  necessary  to  get 
a  Nyquist  plot  of  the  data  (in  rectangular  form) .  The  Nyquist 
plot  is  vector  representation  of  the  real  vs.  imaginary  values 
of  the  data.  On  a  Nyquist  a  pure  mode  would  sweep  out  a 
perfect  circle.  Modes  with  differing  damping  factors  exhibit 
circles  with  varying  radius.  Figure  39  shows  the  Nyquist 
plot  of  the  zoom  data.  Note  there  are  many  modes  present 
with  greatly  varying  damping  factors.  Individual  modes  can 
be  isolated  by  inspection  of  a  small  portion  of  the  rectangu¬ 
lar  data.  For  example.  Figure  40  shows  the  rectangular  plot 
from  about  675  Hz  to  about  850  Hz.  If  the  isolated  portion 
is  converted  to  Nyquist,  as  shown  in  Figure  41,  the  major  mode 
at  772  Hz  is  readily  identified.  The  Nyquist  plots  can  be 
smoothed  (made  more  circular)  by  higher  resolution.  That  is, 
complete  a  zoom  measurement  on  a  smaller  frequency  band. 
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Figure  36.  Zoom  of  test  chamber  section 
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Figure  39.  Zoom  of  test  chamber,  section  1.  Nyquist  plot 


Figure  40.  Zoom  of  test  chamber,  section  1.  Isolated 
portion  of  rectangular  data 


Figure  41.  Zoom  of  test  chamber,  section  1.  Isolated 
portion  of  Nyquist  plot 


All  modes  identified  for  the  test  chamber  and  the  specimen 
(both  fixed  and  free)  are  listed  in  the  tables  included  in 
Appendices  E  and  F. 

It  was  noted  that  there  are  many  more  modes  for  the  test 
chamber  than  for  the  specimen,  but  that  is  expected  because 
the  test  chamber  is  a  much  more  complex  system,  and  the 
boundary  conditions  are  considerably  more  complex. 

The  number  of  modes  for  the  specimen  in  the  test  chamber 
is  found  to  be  26,  whereas  the  number  of  modes  of  the  speci¬ 
men  separated  from  the  test  chamber  is  28.  This  is  good 
agreement,  and  it  should  be  noted  that  the  modes  of  the  test 
chamber  did  not  increase  the  number  of  observed  modes  of  the 
specimen. 

The  average  damping  factor  observed  for  the  test  chamber 
is  about  0.47%,  which  is  a  typical  value  for  plain  carbon 
steel.  The  average  damping  factor  for  the  specimen  (fixed) 
is  about  0.11%,  and  for  the  specimen  (free)  is  about  0.05%. 

As  expected,  the  fixing  of  one  end  of  the  specimen  in  the 
test  chamber  did  increase  the  measured  damping  factor,  but 
with  more  testing  it  is  possible  to  quantify  the  increase 
and  make  suitable  corrections. 

A  plot  of  damping  factors  vs.  frequency  for  the  test 
chamber  is  presented  in  Figure  42.  The  damping  factor  is 
higher  at  lower  frequencies  with  the  maximum  damping  factors 
being  in  the  range  of  1500-3000  Hz.  After  3000  Hz  the  damping 
factor  averages  out  to  about  0.2%  at  higher  frequencies. 
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Figure  42.  Test  chamber — damping  factor  vs.  frequency 


For  the  specimen/  damping  factors  vs.  frequency  (for 
both  fixed  and  free)  are  presented  in  Figures  43  and  44. 

Like  the  test  chamber,  the  damping  factor  appears  to  be 
higher  at  lower  frequency  for  the  fixed  specimen,  and 
without  question  this  same  trend  is  observed  for  the  free 
specimen.  The  fixed  specimen  damping  factor  levels  out  to 
about  0.10%  at  higher  frequencies,  and  for  the  free  specimen 
it  levels  out  to  about  0.025%. 
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Figure  43.  Specimen  damping  factor  vs.  frequency  (specimen  fixed) 


SPECIMEN  (FREE) 


Figure  44.  Specimen  damping  factor  vs.  frequency  (specimen  free) 


VIII.  RECOMMENDATIONS  FOR  FUTURE  WORK 


1.  Expand  the  test  procedure  to  include  random  excita¬ 
tion  of  test  specimen  for  the  energy  input  and  measurement 
of  responses  up  to  15,000  Hz. 

2.  Conduct  similar  tests  on  various  materials  including 
high  damping  materials  such  as  Incramute  and  Sonoston. 

3.  Investigate  the  temperature  dependence  of  viscous 
damping,  especially  for  constrained  layer  and  non-metallic 
composites . 

4.  Investigate  the  changes  in  damping  for  materials 
immersed  in  a  water  medium  at  varying  temperatures. 


APPENDIX  A 


THE  HP-54 51C  FOURIER  ANALYZER 

The  HP-54 51C  Fourier  Analyzer  System  performs  analysis 
of  time  and  frequency  data  containing  frequencies  from  DC 
to  50  KHz.  The  system  analyzes  time-series  data  such  as 
mechanical  vibrations,  sonar  echoes,  tidal  action,  biomedi¬ 
cal  phenomena  such  as  brain  waves  and  nerve  impulses, 
voltages  and  currents  in  electronic  systems,  and  acoustic 
phenomena.  These  analyses  may  detect  signals  hidden  in  noise, 
or  may  locate  critical  functions  in  complex  systems.  Both 
continuous  and  transient  data  may  be  processed. 

Keyboard  programming  allows  the  following  operations 
automatically  without  special  software: 


A. 

Forward  and  inverse  Fourier 

transform 

B. 

Magnitude  and  phase  spectrum 

C. 

Power  and  cross  power  spectrum 

D. 

Transfer  function 

E. 

Coherence  function 

F. 

Convolution 

G. 

Auto  and  cross  correlation 

H. 

Hanning  and  other  weighting 

functions 

I. 

Histogram 

J. 

Scaling 

K. 

Ensemble  averaging  (time  and 

frequency) 
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Six  editing  keys  operate  an  on-line  resident  editor  so 
that  a  sequence  of  steps  configured  into  an  automatic 
measurement  procedure  may  be  changed  on-line  without  the 
need  to  do  off-line  editing,  compiling  and  testing.  In 
fact,  the  series  of  steps  or  programs  used  to  perform  a 
particular  operation  can  be  stored  on  the  Disc  for  easy 
re-entry  into  the  Fourier  Analyzer. 

Data  input  and  output  is  likewise  controlled  from  the 
keyboard.  Data  can  be  entered  in  analog  form  through  the 
four  channel  Analog-to-Digital  Converter,  or  in  digital  form. 

Results  of  all  operations  are  displayed  on  the  oscillo¬ 
scope.  In  addition,  results  can  be  printed  out  in  decimal 
numbers  or  plotted  on  the  Graphics  Terminal  or  an  X-Y  plotter. 

The  Fourier  Analyzer  is  a  completely  calibrated  system; 
all  displays  and  data  output  are  accompanied  by  a  scale 
factor  relating  them  to  physical  units.  This  calibration 
results  from  digital  techniques  being  used  in  all  computations. 


HP-54 51C  SPECIFICATIONS  AND  CHARACTERISTICS 


SPECIFICATIONS 

(Specific  ill  tins  describe  the  rtandard  system's  warranted  performance.) 


ANAIOG-TO-OIGITAL  CONVERTER 

kaput  Voltage  Range:  ±0.12SV  to  ±8V  peak  in  steps  of  2. 

Input  Coupling:  dc  or  ac. 

Input  Channels:  2  channels  wired  lor  4  standard.  4  channels  op¬ 
tional  with  piug-in  cards. 

Resolution:  12  bits  including  sign. 

Input  Frequency  Range:  dc  to  SO  kHz,  5  Hz  to  50  kHz,  ac 
coupled  (100  kHz  optional!. 

Simple  Rale: 

Internal:  100  kHz  max.  d.  2.  3.  or  a  channels  simultaneously!. 
(200  kHz  optional  on  1,  2,  3.  or  4  channels,  i 
|S0  kHz  max.  (3  or  4  channels  simultaneously!. t| 
External:  An  external  time  base  may  be  used  to  allow  external 
control  of  the  sampling  rate  up  to  100  kHz  1200  kHz 
optional).  One  sample  can  be  taken  for  each  clock 
pulse  <TTl  level). 

Internal  dock  Accuracy:  ±0.01%. 

OtSPlAY  UNIT 

Vertical  Scale  Calibration:  Data  in  memory  is  automatically  scaled 
to  give  a  maximum  on-screen  calibrated  display.  The  scale  factor 
is  given  in  volts/division,  voltst/divition.  or  in  dB  offset. 

Leg  Display  Range:  80  dB  with  a  scale  factor  ranging  from  0  to 
+99R  dB.  Offset  selectable  in  4  dB  steps. 

linear  Display  Range:  ±4  divisions  with  scale  factor  ranging 
from  1  X  10-t>2  to  S  x  lOStz  in  steps  of  1.  2.  and  5. 

Digital  UP/DOWN  Scale:  Allows  8  up-scale  and  2  down-scale 
steps  (calibrated  continuous  scale  factor! 

Horizontal  Scale  Calibration: 
linear  Sweep  length:  10, 10.24  or  12.8  divisions, 
log  Horizontal:  0.5  decades/division. 

Mathew:  Intensity  markers  every  8th  or  every  32nd  point. 

RASE  SOFTWARE 

Ttawdetm  Accuracy:  The  expected  rms  value  of  computational 
error  introduced  in  either  the  forward  or  inverse  FFT  will  not 
exceed  0.1%  of  the  rms  value  of  the  transform  result. 

Dynamic  Rangt:  >  75  dB  for  a  minimum  detectable  spectral  com¬ 
ponent  in  the  presence  of  one  full  scale  spectral  component 
after  twenty  ensemble  averages  for  a  block  size  of  1024. 


EXECUTION  TIMES* 

Fourier  Transform:  <55  ms 

Stable  Power  Spectrum  Average:  <80  ms 

Stable  Tri-Spectrum  Average:  <  220  ms 

REAL  TIME  BANOWIDTHS* 

Fourier  Transform:  >7.5  kHz 

Stable  Power  Spectrum  Average:  5.4  kHz 

Stable  Tri-Spectrum  Average:  1.9  kHz 

MASS  STORAGE  SOFTWARE 
MAXIMUM  REAL  TIME  DATA  ACQUISITION  RATE 
(Single  Channel): 

BS  25*:  10  kHz 

BS  1024:  39  kHz  (25  kHzt) 

BS  409*:  80  kHz  1 30  kHzt) 

OFF-LINE  BSFA  SOFTWARE 

Center  Frequency  Range:  dc  to  one-half  the  Real  Time  Data 
Acquisition  Rate. 

Center  Frequency  Resolution:  Continuous  resolution  to  the  limit 
of  the  frequency  accuracy  for  center  frequencies  >0.02%  of  the 
sampling  frequency. 

Frequency  Accuracy:  ±0.01% 

Bandwidth  Selection:  In  steps  ol  f/Sn  where  n  »  2,  3.  4.  etc. 

Max.  Resolution  Enhancement:  >400 

Dynamic  Range:**  90  dB  from  peak  out-of-band  spectral  com¬ 
ponent  to  the  peak  level  of  the  passband  noise. 

B0  dB  from  peak  in-band  spectral  component  to  the  peak  level 
of  the  passband  noise. 

Out-of-Band  Refection:  >90  dB 

Famband  Flatness  of  the  Digital  Fiber:  *0.01  dB 

ENVIRONMENTAL  CONDITIONS 
Temperature  Range:  0*C  to  40*C  <104*F|. 


•*or  bond  hmOOd  roMtom  new  typo  NgEM*  m  bfac*  Wf  1024.  no  dupiov.  no  Honntng 
orght  «»EH|ii  o i  t  po*nf  spacmufi  «  Mocft  ur*  tOJ4  Hockocod  by 

TO  d*  M  n*ci  com or  of  rto  bond 

'.Thom  earn  WO*y  to  lyttwra  ~*h  moduht  S4MI  ond  $4451 A/I  htwng  4  tonol  profit 
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SUPPLEMENTAL  CHARACTERISTICS 
(Supplemental  Characteristkl  are  intended  Is  proride  useful  miormudon  for 
system  appttcadoa»  by  giving  lypkai  but  not  warranted,  periormance  parameters.) 


ANALOG-TO-OtCITAl  CONVERTER 
Input  Imped tncei  1  MO  in  parallel  with  <75  pi. 

lampli  Rate  Conduit 

Maximum  frequency  Mode:  Maximum  frequency  selectable, 
from  0.1  Hz  lo  50  kHz  (100  kHz  optional!  in  steps  of  1,  2.5,  5. 
This  mode  automatically  sets  maximum  frequency  indepen¬ 
dent  of  block  size. 

frequency  Resolution  Mode:  Frequency  resolution  selectable 
from  0.2  mHz  to  1000  Hz  in  steps  of  1.  2,  5.  This  mode  auto¬ 
matically  sets  frequency  resolution  and  sample  record  length 
independent  of  block  size. 

Input  Model  There  is  a  buffered  and  non-buffered  analog  mode. 
In  the  buffered  mode,  other  operations  can  be  performed  on 
previously  collected  data  while  the  ADC  collects  current  input 
data  into  a  buffer. 

DISPLAY  UNIT 

Data  may  be  displayed  in  single  sweeps  or  refreshed  continu¬ 
ously  in  the  following  forms: 


Amplitude 
Real  Pan 
Imaginary  Part 
Magnitude  i  Unear  or  Log) 
Phase 


Time  (Linear  or  Log) 
frequency  (Linear  or  log) 
frequency  (Linear  or  Logi 
Frequency  (Linear  or  Log) 
frequency  (Linear  or  Log) 


j  Imaginary  Pwt  (NVQUIST  PLOT)  Real  Pan 

CRT  Positioning:  Three  markers  to  aid  in  adjusting  trace  posi¬ 
tion  as  well  as  venical  and  Horizontal  controls  are  provided  for 
display  positioning. 

Origin:  Left  edge  of  display,  zero  amplitude. 

♦fit  Positive  full  scale,  center  of  display. 

-fS:  Negative  full  scale,  center  of  display. 

Analog  Hotter  Output:  Any  displayed  data  can  be  output  to  a 
plotter  or  remote  oscilloscope.  HP  106408  interface  is  required 
for  operation. 

Ampttude:  0.5V  per  oscilloscope  display  division. 

Uneasily:  0.1*  full  scale. 

Interpolation:  Linear  interpolation  in  0.01%  steps. 

Type  of  Display:  Points,  bars,  or  continuous  (interpolation:. 

RASE  SOFTWARE 

System  Accuracy  and  Range:  The  Fourier  Transform  is  imple¬ 
mented  using  conditional  scaling  for  maaimum  accuracy  with 
no  data  overflows  allowed.  All  calculations  use  floating  point 
arithmetic  on  a  block  basis  with  full  16-  and  32-bit  arithmetic 
where  applicable. 

Data  Ward  Size:  16-bit  imaginary  with  32  bits  preserved  for 
double  precision  functions.  Division,  addition,  or  subtraction 
operations  performed  in  16  or  32  bits  depending  on  data. 
Maximum  Block  Size:  4096  time  domain  points. 

Minimum  Slock  Size:  M  time  domain  points. 

Penasansnt  Data  Space:  28K  words  >  16K  words  standard  in  systems 
with  serial  prefix  below  1842.  optionally  expandable  to  28K  with 

option  0111. 

Nrmenent  Program  Space:  32K  words. 

RSFA  SOFTWARE 

Maeknum  BSFA  Hocksize:  1024  time  domain  points.  1 2048  with 
option  670.)  See  also  Table  5-1  in  BSFA  section  of  manual. 


MASS  STORAGE 

Diac  Unit: 

Capacity:  2.45  megawords 

Data  Transfer:  2.5  million  bits/second 

Discs:  2  |1  fixed,  1  removable) 

%  of  Real  Time  at  16t  kHz  ADC  Sampling  Rale  (Single  Channel): 

IS  256:  10%  RS  1024:  39%  (25%fl  BS  4096:  80%  i30%ti 
Number  of  Records  Per  F0c: 

Data  Slock:  214  records  >4096  blocksize  maximum/recordl. 
ADC  Throughput:  199  records  .40%  blocksize  max./record  i. 
Program  Stack:  138  records  i470  steps/recordi. 

ASCII  Test:  690  records  1 128  words/record'. 

Index:  69  records  i10  pointers/ record  I. 

System  Cordoad:  4  records  (32K  words/recordi. 

Common:  286  records  (2S6  words/recordi. 

Overlay:  20  overlays  SK  words  maximum,  7K  words  maximum 
with  option  261  or  265). 

OPTIONS  261  R  26S 

The  magnetic  tape  options  are  used  for  ADC  Throughput  only. 

Maximum  Real  Time  Data  Acquisition  Rate  (Single  Channel): 


Opt  261  Opt  265 

BS  2S6:  6  kHz  BS  256:  9  kHz 

BS  1824:  12  kHz  BS  M24:  21  kHz 

BS  4096:  15  kHz  RS  4096:  30  kHz 

klumher  of  Tracks:  9 
Read/Write  Speed:  4S  ips 

Density.*  Option  261.  800  bpi.  Option  265,  1600  bpi. 

Data  Transfer  Rate:  Option  261.  36K  cps  max..  Option  26S.  72K 
q»  max. 

Rack  Height:  610  mm  i24  in.i 


POWER  REQUIREMENTS.  SIZE.  WEIGHT 
Power  Source:  115/230  volts  ±10%.  50/60  Hz.  1800  wans  tvpical 
for  base  system. 

Size:  Dimensions  are  for  a  typical  system  (excluding  cabinet  and 
rerminall. 

Height:  771  mm  i28  m.i;  Width:  425  mm  <16Va  m.i; 
Depth:  616  mm  i24Va  in.). 

Cabinet: 

Panel  Height:  1422  mm  (S6  in.i 

Overal:  1631  mm  64</a  in. i.  height;  533  mm  i21  in. i,  width; 
762  mm  (30  in.i.  depth 

Weight:  Net  weight  for  a  base  system  (excluding  terminal  (  — 
163.3  kg  (358  Ibs.i. 

Price  and  Ordering  Information:  Consult  the  5451C  Fourier 
Analyzer  System's  Ordering  Information  Guide. 

SYSTEM  INSTALLATION 

Included  in  the  5451C  System  is  on-site  installation.  On  installation, 
a  trained  Hewlett-Packard  representative  will  perform  an  opera¬ 
tional  demonstration  to  ensure  that  the  system  is  functioning  normally 

TRAINING 

A  course  on  Fourier  analysis  and  system  operation  is  optionally 
available  at  HP's  Santa  Clara.  California  facility.  On-site  training 
can  also  be  provided,  if  desired. 

i  thru  awxxwifw  wpty  (o  W*h  iwoZuiaa  S4SSZ  an*  saastA/a  Aawwj  a  ixnal 

prxtw  lower  than  taaz 
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APPENDIX  B 


HP-5451C  TRANSFER  FUNCTION  AND  POWER  SPECTRUM 
FLOW  CHART  AND  PROGRAM  LISTING 


This  appendix  contains  flow  charts  and  listings  of  the  pre-written  soft  key  programs  shipped  on  the  5451C 
operating  disc  pack.  These  programs  are  for  the  Gold  Key  F2  (Transfer  Function)  and  F5 1 Power  Spectrum) 
soft  keys.  The  purposes  for  including  these  flow  charts  and  listings  are  as  follows: 

1.  For  those  who  wish  to  modify  the  programs  to  better  fit  specific  applications. 

2.  As  an  example  for  those  who  wish  to  write  their  own  soft  key  programs. 

3.  As  a  model  for  re-entering  portions  that  may  have  been  accidentally  written  over  on  the  disc 
(programs  reside  in  unprotected  areas). 

As  an  aid  toward  rapid  understanding  of  these  programs,  a  certain  programming  "style”  has  been  used.  As 
a  result  of  this  style,  the  programs  are  longer  than  they  would  otherwise  need  to  be,  but  are  easier  to 
comprehend.  Some  elements  of  the  style  are  described  below. 


BLOCK  STRUCTURE  THROUGH  USE  OF  LABELS 

The  code  in  the  programs  is  organized  in  “blocks"  —  functional  segments  which  are  delimited  by  LABEL 
instructions.  A  convention  is  followed  in  the  choice  of  label  numbers.  The  beginning  of  each  block  is 
designated  by  a  label  number  ending  in  a  multiple  of  50,  e.g.  L1050.  The  block  ends  with  the  label  9  higher, 
e.g.  L1059.  Within  the  block,  label  numbers  are  between  these  limits,  e.g.  L105I,  L1052  etc.  One  block  may 
contain  others  —  each  will  use  the  above  delimiting  convention. 

Although  the  use  of  these  labels  lengthen  the  programs,  it  makes  them  much  easier  to  understand,  and  to 
correlate  with  the  flow  diagrams.  For  example,  when  several  branches  of  a  program  go  to  the  same  point, 
each  branch  will  jump  to  its  own  appropriate  label.  Since  these  labels  all  appear  together  at  one  point  in 
the  program,  one  can  tell  that  multiple  branches  converge  at  that  point.  An  example  is  in  the  Transfer 
Function  program,  stack  54,  where  labels  1209, 1259  and  1609  all  appear  together.  This  is  just  one  example 
of  the  way  in  which  complex  internal  linkages  in  the  program  are  made  more  visible. 


BRANCHING  THROUGH  USE  OF  “COMPUTED  GOTOV’ 

In  most  complex  programs,  branching  is  common.  One  means  of  branching  is  to  use  an  IF  statement, 
provided  in  Keyboard  language  by  the  "GOLD  KEY"  "SKIP"  instruction.  When  there  are  more  than  two 
possible  branches,  however,  use  of  IF  branching  tends  to  get  complicated,  involving  multiple  decision 
points.  The  "Computed  GoTo”  or  "switch"  type  of  branching  statement  is  more  suitable  in  such  cases  for 
simplicity  of  understanding.  It  has  been  used  extensively  in  these  keyboard  programs,  even  for  simple  two 
way  branching.  By  standardizing  on  it,  the  code  becomes  recognizable  and  easier  to  read. 
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In  Keyboard  language,  the  “Computed  GoTo”  is  implemented  by  computing  the  number  of  a  label  and 
jumping  to  it.  The  following  is  an  example. 


l  lisa 

Starr  of  branch  block 

Y  A*  1  1ISI  11D 

Sat  variable  0  to  1050 
♦  (the  valoe  of  variable 
It) 

J  ID 

Imp  to  the  label  nmber 
in  variable  1 

L  1IS1 

# 

« 

Code  to  be  executed  if 
variable  11  =  1 

J  1059 

Go  to  end  of  block 

L  10S2 

Code  to  be  executed  if 
variable  11  *  2 

L  IDS?  End  of  branch  black 


USE  OF  SUBROUTINES 

In  these  programs  there  are  several  functions  which  have  been  set  up  as  subroutines.  These  include  the 
parameter  entry  routine  ( see  below)  and  the  measurement  routines.  The  measurement  routines  are 
handled  this  way  to  simplify  the  flow  of  the  main  program  and  allow  easy  replacement  of  the 
measurement  code  in  case  you  modify  it  in  some  way. 


PARAMETER  ENTRY  ROUTINE 

Since  these  programs  ask  the  operator  for  many  input  parameters,  a  single  subroutine,  LABEL  100  in 
STACK  0,  handles  all  the  parameter  entries.  The  routine  is  called  with  variable  parameters  1  and  2  equal  to 
the  lower  and  upper  limits  on  the  range  of  allowable  operator  inputs.  The  routine  reads  your  input, 
checks  it  against  the  range  limits,  and,  if  it  is  valid,  passes  it  back  to  the  calling  program  in  (floating  point) 
variable  parameter  2000.  If  the  input  is  out  of  range,  this  routine  notifies  you  and  waits  for  the  new  input. 
The  routine  will  not  return  to  the  cajiing  program  until  a  valid  input  has  been  received. 


PRECAUTIONARY  NOTES 

The  following  precautions  apply  to  the  operation  of  the  preprogrammed  measurements: 

1.  When  using  the  standard  software  zoom  fBSFA),  the  measurement  blocksize  can  be  no  larger  than 
1024.  When  using  the  Option  670  Fourier  Preprocessor  for  BSFA  measurements,  the  maximum 
blocksize  is  2048. 

2.  The  maximum  center  frequency  you  may  specify  for  a  BSFA  measurement  is  32767  Hz. 

3.  The  messages  CF  WHAT?  or  BW  WHAT?  may  result  if  the  center  frequency  and/or  bandwidth  you 
have  chosen  for  your  measurement  are  such  that  the  BSFA  analysis  band  is  either  less  than  0  or 
greater  than  the  AOC  Fmax  setting.  Specifying  different  parameters  should  remove  this  problem. 
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4.  The  message  DL  WHAT?  may  occur  when  performing  the  on-line  BSFA  measurment.  This  is 
because  the  display  is  active  during  the  on-line  measurement.  To  remove  this  problem,  either 
reduce  the  measurement  bandwidth  (thereby  increasing  the  zoom  power  and  lowering  the  data 
rate  into  the  computer),  or  edit  the  appropriate  keyboard  stacks  (stack  56  for  transfer  function, 
stack  61  for  power  spectrum)  to  remove  parameter  n3  from  the  calls  to  User  Prog  45  for  the  on-iine 
measurement  (refer  to  commented  program  listings  which  follow). 

5.  When  performing  an  off-line  BSFA  measurement  with  an  optional  Mag  Tape  unit,  perform  the 
following  steps  before  making  the  measurement. 

Set: 

ADC  SAMPLE  MODE  to  INTERNAL  KHz//*s 
MULTIPLIER  to  100/10/5 
INPUT  SELECTOR  to  A 
TRIGGERING  to  FREE  RUN 
OVERLOAD  VOLTAGE  A  to  CHECK 

Enter: 

BLOCKSIZE  4096  ENTER 
MASS  STORE  32  ENTER 

MASS  STORE  22  SPACE  1  SPACE  ISO  ENTER 
MASS  STORE  32  ENTER 

This  writes  150  records  of  data  on  the  magtape  so  that  the  magtape  will  be  able  to  position  to 
record  135  on  the  tape  when  the  ADC  throughput  is  performed.  It  will  position  by  looking  for  the 
interrecord  gaps  written  by  the  WRITE  ADC  throughput  command. 

6.  After  completing  a  BSFA  measurement,  be  sure  that  all  data  space  declared  by  the  zoom  programs 
is  released  by  pressing  RESTART. 

As  you  go  through  the  flow  charts  and  commented  listings,  remember  that  these  are  only  examples  of 
programming  the  soft  keys  FI  through  F6  on  the  Keyboard.  It  is  up  to  you  to  determine  which,  if  any, 
portions  of  these  programs  should  be  maintained.  Because  these  programs  are  stored  in  unprotected 
areas  of  the  Disc,  there  is  the  possibility  they  can  be  written  over.  If  this  should  happen,  you  should  enter 
the  program  stacks  from  the  listing,  substitute  your  own  program,  or  copy  from  your  back-up  disc. 

The  soft  key  programs  and  the  associated  ASCII  text  and  variable  parameters  were  originally  stored  on  the 
system  disc  pack  in  Files  3,4,  and  7.  The  records  used  are  as  follows: 

File  3  .'Keyboard  Programs)  File  4  (Text  Buffersi 

Record  0  Text  buffers  51  through  55 

Records  51  through  62  ASCII  records  32‘f5"through3‘W* 

File  7  (Common) 

Common)  Record  0 

•This  assumes  thar  there  areJTOO  records  in  File  4.  If  not,  the  first  and  last  ASCII  text  records  may  be 
computed  as  follows: 

First  record  number  =  NR  -  (5  x  last  text  buffer  number) 

Last  record  number  =  NR  -  (5  x  first  text  buffer  number)  +  4 
where  NR  =  number  of  records  in  File  4. 


The  allocations  above  should  be  kept  in  mind  so  the  above  records  are  not  inadvertently  altered  or 
destroyed  when  using  the  Fourier  system.  Should  you  wish  to  alter  the  allocations,  you  will  also  have  to 
modify  the  keyboard  programs  to  reflect  such  changes. 


ii 
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TRANSFER  FUNCTION 
FLOW  DIAGRAM 


H 


LABEL  2000 


RANDOM - 


SET  # AVERAGE 
♦  I  LI  2001 


^  TYPE  ^ 
.MEASURE 


TRANSFER 

NyFUNt/ 


TURN  ON  HF  DAC 


•  TRANSIENT 


'BASEBAND 


SELECT  TYPE  BSFA 
AND  CF/8W 


NUMERICAL 


\  ..  -  ■<  CFIBW  - -  j 

Vi  (LI 301  rV^TRY^ (LI 302) 


CURSOR 


BASEBANO  TRANSFER 
FUNCTION 
MEASUREMENT 


ENTER  CF/SW 
FROM  TERMINAL 


SET  BW 
WITH  CURSOR 


WAS  ^ 
DAC  ON 
•v  f  / 


OLO  ^ 
(L1461) 


TURN  OFF  DAC 


1  (LI  551) 

[abort > 


USE  LISTEO 
TRACKS 


SET  THRU-PUT 
TRACKS 


off-line  , 

SOFTWARE 


,  ON-LINE  PRE-PROCESSOR 


TYPE  \ _ 

BSFA 

N.  OFF-LINE 

PRE-PROCESSOR 


_ I— — — 

INITIALIZE  BSFA 

INITIALIZE  BSFA 

INITIALIZE  BSFA 

OFF-LINE 
SOFTWARE  BSFA 

OFF-LINE 

544 70A  BSFA 

ON-LINE 

54470A  BSFA 

- - - - 

(L1309.  L1269.  LlSOBI 


DISPLAY  BLOCKS  S 
PLOT  IF  DESIRED 


PRESS  CONTINUE 


CONTINUE 
V  f  X 
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L 


0 


STACK  S 

Sift  Hi  juno  dirtc*arv 
■GOLD  \EV*  ‘Fn*  causes  juna  to  LABEL  i 
in  this  stock.  Use.'  can  aragran  tit  staci 
anywhere  fnn  'here. 

Transfer  Function  ms  Piuer  Spectrun 
art  arepregranned  in.  Tntse  nay  also 
bt  cncngec, 


L 

2 

F2  TRANSFER  FUNCTION 

j 

1800 

SI 

Label  1300  in  stack  51 

L 

5 

F5  POkfER  SFECTRJh 

J 

1000 

57 

Label  1030  .rt  stack  57 

FI  USER  ASSIGNABLE 

L 

1 

L 

3 

F3  USER  ASSIGNABLE 

L 

4 

F4  USER  ASSIGNABLE 

L 

6 

F6  USES  ASSIGNABLE 

D 

Pause  here  if  «r..:ss. cried  fur.cuin  called 

J 

6 

Loop  if  ‘CONTINUE*  pressec 

1. 

100 

PARAMETER  ENTRY  RGUTINE 

Y  t 

2300 

Rtad  user  entry  *'.t 3  .if. ail*  2033 

r  if 

2003 

ID  S  1 

Check  :**  entry  '*  it-er  i.n.t 

Tea  snail.  Cue  .jerattr  w»tn  *anqe 

Y  J 

99 

1 

Y  P 

i 

Y  * 

96 

\ 

Y  P 

■j 

•» 

; 

136 

*1 

and  loop  for  new  entry 

L 

ill 

Y  I? 

2000 

2D  5  -1 

Check  >'  ent-y  •.*  jpper  Lju' 

r  i> 

09 

1 

Tio  iargt.  lie  aperctar  with  -ar.ce 

Y  P 

* 

Y  U 

98 

1 

Y  ? 

2 

J 

166 

m* 

and  loop  e ^  new  e  .**v 

( 

Return  to  coi-.r:;  ?;:rcn“-entry  ;n 
uariatle  2010*  wi'tin  sjecifiea  rarue 

L 

189 

:s: 
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"JfL 


L 

-51 

L 

1140 

MS 

37 

MS 

17 

T 

5338 

51 

Y  W 

1 

1 

Y  U 

2 

1 

Y 

1 

l 

Y  ‘ 

2 

3 

J  " 

100 

0 

Y 

11 

2000D 

L  ‘ 

1050 

Y  At 

9 

1050 

J 

9D 

L 

105i 

L 

1053 

J 

105? 

L 

1052 

Y  4 

3 

i 

Y 

i 

3 

Y  ' 

2 

9999 

J  ’ 

100 

9 

Y 

12 

20I0D 

t  * 

U 

1059 

Y  4 

5 

1 

BS 

512 

D 

Y  U 

& 

1 

Y 

i 

9 

Y  ‘ 

2 

i 

J  ’ 

199 

0 

Y 

14 

2000D 

L  ’ 

1100 

Y  A+ 

Q 

1101 

J 

OD 

L 

1101 

J 

1199 

1 

W 

1102 

Y 

100 

i) 

L 

1109 

Y  4 

n 

7 

i 

D 

Y 

5838 

52 

L 

1150 

Y  A+  0 

1150 

J 

0D 

L 

1151 

Y  W 
J 

8 

1 

1159 

L 

1152 

B 

6 

12C 

B 

6 

J 

US? 

L 

1153 

Y  U 

9 

l 

L 

115? 

1 

Y  « 

10 

RA 

F 

TL 

I 

2100 

1 

STACK  51 

TRANSFER  .-UNCTION  PROGRAM 
Recall  previous  variable  parameters 

Read  text  buffer  951 
Print  nessages  i  4  2 

Set  allowable  range  'or  user  entry 

Call  user  entry  routine  LUO,  stack  3 


BLOCK  1050 

liD  Branch  per  value  of  vgr,able  ii 
to  LiuSl,  1052  or  1053 


Parameter  entry 
Range:  I  *  1 
Variable:  l2 

END  OF  BLOCK  1050 

Block  sire  512— Jser  ,iay  change  after  pause 
Pause  for  operator  action 

Paraneter  entry 
Range:  0  -  t 
Variable:  14 

BLOCK  1100 
14D 


-2  Set  2  filters— auto  select  node 
END  OF  BLOCK  1100 

Pause  for  operator  action 

BLOCK  1150 

HD  Branch  per  variable  :: 
to  LllSi,  1152  or  US3 


Fill  DAC  buffer  and 
turn  in  DAC 


END  OF  BLOCK  1150 

Analug  input  in  ’REPEAT*  tu  nonitor  inputs 

Fewer  spectrun  (log)  for  user  reference 
1  Junp  to  next  stack  and  continue 
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L 

-52 

STACK  52 

L 

2118 

LOO?  POINT  r0R  TRANSFER  "UNCTION 

Y 

5838 

52 

Y  tf 

11 

1 

Y 

1 

1 

Poroiieter  entry 

Y  “ 

2 

32787 

Range:  1  -  32787 

J  * 

181 

0 

Variable:  15 

Y 

IS 

20000 

L  ‘ 

i230 

BLOCK  1200  (enos  .n  track  54) 

Y  A* 

8 

1200 

110 

Branch  per  vcriohle  li 

J 

80 

tt  1.1201,  12C2  or  1203 

L 

1282 

* 

y  y 

22 

1 

3 

Paese  for  operator  action 

B 

8 

Torn  on  DAC 

J 

3808 

3 

1 

Svbrovt.ne  call  -  baseband  neasurenent 

B 

T»rr.  off  DAC 

J 

120? 

2 

1 

then  30  to  er.a  of  block  :230 

L 

1203 

y  u 

55 

4 

A 

0 

Pause  for  operator  act. on 

w 

3030 

3 

4 

Sobrootine  call  -  baseband  reasjrenent 

T 

1# 

1209 

2 

; 

then  40  to  erd  of  slock  1230 

L 

1201 

y  y 

12 

1 

Y 

1 

i 

Paraneter  entry 

Y  * 

2 

2 

Range:  1  -  2 

J 

100 

0 

Variable:  18 

Y 

18 

2500D 

L  ’ 

1250 

BLOCK  1250  (ends  in  stack  54) 

Y  A+ 

3 

1250 

160 

Branch  per  .'enable  it 

J 

CD 

to  L12S1  or  1252 

L 

1251 

Y  U 

22 

\ 

3 

Pause  for  operator  action 

J 

3000 

3 

i 

Subroutine  call  -  baseband  neasorerer 

r 

y 

1259 

2 

1 

then  go  to  end  of  sloe*  1250 

L 

1252 

Y  U 

13 

i 

•/ 

f 

i 

» 

i 

Paraneter  entry 

Y  ' 

2 

3 

Range:  1-2 

J 

100  “ 

0 

Variable:  17 

Y 

17 

203QC 

Y  _ 

5336 

53 

Y  U 

IS 

i 

i 

1 

Parameter  enTry 

v  ** 

2 

£ 

Range:  1  -  2' 

t  “ 

100 

0 

var.abie  :3 

/ 

* 

IS 

290GD 

L 

1300 

3L2CK  1300  (ercs  ..a  :*ack  32) 

;  A* 

0 

1300 

:8D 

Branch  oer  var.asle  i2 

J 

0D 

to  L:301  or  i3C2 

w 

1301 

y  y 

18 

4 

V 

1 

4 

Parameter  er.-v 

Y  * 

2 

32787 

Range:  i  -  32787 

10G 

i# 

Variable:  i? 

Y 

1? 

200 'it 

Y  S 

17 

e 

Y 

1 

i 

Parameter  entry 

Y  3- 

5 

fa 

l°0 

1 

Since.  1  -  {value  of  oar  19  -  c. 

* 

y 

100 

0 

Variable:  20 

V 

i 

20 

2CC0D 

4 

1309 

1 

* 

L 

1302 

J 

1303 

i 

1 
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L  -S3  STACK  S3 

L  1313 


Y  U 

18 

1 

Cirser  an.  Cursor  paraneters  to 

/. 

2111 

variables  2000  -  2002  (200isfrequency) 

Y 

2103 

20010 

Variable  2003  gets  ist  cursor  frequency 

Y  0 

1? 

1 

/. 

2101 

Cursir  in,  paraneters  to  2000  -  2002 
Cursor  #ff 

/. 

-1 

Y 

2004 

21019 

Variable  2004  qets  2nd  carsir  frequency 

L  " 

1350 

2003 

block  i350 

Y  IF 

20049 

2  -2 

If  let  cursor  frequency  >  2nd 

Y 

Y  ' 

2004 

2103 

20039 

20019 

then  swap  then 

L  “ 

13S9 

EN9  OF  BLOCK  1350 

Y  A- 

2030 

20G4D 

20039 

Y 

Y  T 

20 

2000 

29009 

Variable  20  qets  :«*n  bandwidth 

Y  A* 

2001 

20039 

20009 

Y 

L  “ 

19 

1410 

20119 

Variable  i9  sets  zoon  center  frequency 
BLOCK  1400 

Y  IF 

21 

199 

1  -2 

If  bandwidth  not  <  ctr  treq 

Y  A- 

20 

199 

1 

then  bandwidth  =  err  freq  -  1 

EN9  OF  BLOCK  1400 

L 

1409 

L 

1309 

END  OF  BLOCK  1300  vfron  stack  52) 

L 

14S0 

BLOCK  1450  (ends  in  stack  54) 

Y  At 

0 

1450 

170 

Branch  per  variable  i7 

J 

L 

0D 

ti  L1451,  1452  ur  1453 

14S1 

Y 

5838 

54 

Y  J 
! 

22 

1459 

1 

1 

1 

Pause  fir  operator  action 

1452 

L 

14S3 

y  v 

20 

1 

Y 

1 

1 

Paraneter  entry 

V  - 

2 

2 

Range:  1-2 

X  ' 

100 

a 

Variable:  24 

Y  . 

24 

20009 

ft* 

1500 

BLOCK  1500  (ends  in  stack  54) 

Y 

5838 

54 

r  At 

a 

1500 

240 

Branch  per  variable  24 
to  L1S01  or  1502 

i 

ID 

L 

1501 

J 

1509 

♦ 

1 

L 

1S12 

y  y 

28 

1 

Y 

1 

1 

Paraneter  entry 

Y  " 

2 

3 

Range:  1-3 

;  • 

100 

0 

Variable:  2000  (tenporarvi 

l 

1551 

BLOCK  1550  (ends  in  stack  54> 

Y  At 

0 

1550  ’ 

20000 

Brancn  per  variable  2000 
to  L1551,  1552  or  1553 

J 

09 

L 

1551 

D 

lSSi 

Pause  b  loop  to  here-abort 

I 

1 

ft* 

1552 

J 

1554 

1 

1 

L 

1553 

J 

1555 

t 

* 

1 

t 


I 
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L 

-54 

L 

1554 

Y 

22 

134 

Y  ‘ 

21 

64 

*  m 

4 

1559 

L 

1555 

Y  N 

27 

1 

Y 

i 

1 

Y  ‘ 

2 

197 

r  ■ 

109 

0 

Y 

22 

2090D 

Y  S- 

21 

198 

■ 

1559 

NS 

32 

22D 

Y  11 

22 

4 

* 

D 

Y  BS 

13 

BS 

2048 

NS 

22 

2 

BS 

13D 

Y  U 

21 

1 

L 

150? 

L 

L 

Y  A+ 

0 

1696 

STACK  54 


Var.oble  22  gets  default  start  track 
loanable  21  gets  default  9  of  records 


Parameter  entrv 
Range:  1  -  197 
Variable:  22 

Variable  22  aets  user's  start  'rack 
variable  21  bets  ♦  a*  records  left 
END  CF  BLOCK  1550  1?>m  stack  3) 
Position  throughput  file  tc  start  tract 


Pause  for  operator  action 
Variable  13  gets  curren»  block 
Set  blec«  site  to  *ar  for  taro 


Perform  AX  throughput 
Restore  slock  size 


ock  Size 
.troughput 


END  OF  SLOCK  1505  ,fro«  srac*  53) 
END  OF  BLOCK  1453  (fron  stack  S3) 
BLOCK  1650 

Branch  per  variable  :? 


; 

JD 

ti  LlfeOl ,  1602  ur  io03 

L 

1601 

Y 

40 

19D 

20D 

Ini'ializt  zoon—on  li.te  preprocessor 

J 

4900 

2 

[ 

Subroutine  call-zoom  measurement 

Y 

40 

3 

Reset  zoom  to  baseband 

J 

1609 

Ge  to  end  of  block  1601 

L 

1602 

Y 

43 

19D 

29D 

22D 

InitiaLze  zuoit—cff  line  preprocessor 

J 

4500 

2 

1 

Subroutine  call— zoom  measurement 

V 

j 

43 

0 

Reset  zoon  to  basebanc 

1609 

Go  to  end  of  block  16(1 

L 

1603 

Y 

41 

19D 

200 

22D 

Initialize  zoom— off  lire  software 

J 

4500 

*3 

Subroutine  call— zoom  measurtmen* 

V 

41 

0 

Reset  zoom  to  dasesar.d 

L 

1609 

END  OF  BLOCK  1690  stack  uD' 

L 

1259 

END  OF  ?LGCX  125?  ■.from  stac<  52) 

L 

120? 

END  OF  BLOCK  1209  ifron  stack  52/ 

Y 

5838 

54 

Transfer  function  c^nplete 
Print  .nstructions  to  operator 


Save  variable  parameters 

Pause  for  operator  action 
1  .  Loop  to  repeat  transfer  function 

END  OF  TRANSFER  FUNCTION  teicept  susroutmes) 


I  -ss 
l  3m 
Y  5838 
t  b 


I  3(18 

a  2 


CTATY  CC 

BASEBAND  TRANSFER  FUNCTION  SUBROUTINE 

Check  ADC  Mput  selector.  If  net  2  channel 
infer*  operator 

savae  fnr  operator  to  ctrrect  it  ana 
loop  to  check  again. 

Clear  blacks  needed  for  averaging 


}  3858 


i 

3859 

L 

3IS3 

Y  H 

29 

•» 

3854 

Y  i 

31 

BA 

8 

Y 

5 

F 

« 

a 

8 

a 

l 

SP 

8 

* 

33S4 

i 

9858 

CH 

8 

l 

3159 

TL 

0 

Tl 

2 

BLOC!  3858  ,  * 

11D  Branch  per  'liable  11 
to  .3851.  3852  er  305.' 

L3851  h  3852  eg#i»aler.x— non-transient 

Analog  input 
Nairn 

Transfer* 

Clear  DC 

2  Average  spectra 

8  Loop  for  specified  *  :f  averages 

\  Subroutine  call  -  clear  2nd  nclf  bioc» 

2  Cs*oote  transfer  function 

Correct  anplitode  for  Hanning 


l3853  transient  signals 
Instroct  eperater  to  .-span  on  cse 

Coe  with  'beep* 

Analog  .np«x 

Check  for  overload— it  so,  coe  again 
Transfer* 

Clear  DC 

Aoerace  spectra 

Loop  for  specified  8  .t  overages 
Subroutine  call  -  ciear  2nc  holt  bleu/ 
Connote  transfer  function 
EMO  OF  BLOW  al-50 
Take  logs  of  spectra 


Return  to  calling  procrc* 


L  -54 

L  4888 

Y  5838 

Y  4 

31 


J 

a 

a. 

CL 

a 

L 

Y 

S t 
8 

CM 

* 

4 

t 

L 

X< 

t) 

X( 

X/ 

x< 

X) 

x< 

x> 

X( 

ft 

TL 

TL 

< 

L 

a 

CL 

a 

CL 

w 

Y 

Y 
8 
SP 

# 

CM 

J 

< 


4181 

7 

8 
9 
18 

4118 

45 

5 

4818 

5 

4288 


4288 

9 

4 

8 

3 

7 
2 

4 
1 

5 

8 
2 

3 

4 


4588 

7 

8 

9 

10 
4518 

45 

4! 


55 

2 


S 

4518 

S 

4218 


iSD 

2 


4 

5 


4 

2 

150 

2 


SUBROUTINE  FOR  ON  LINE  TOON  TRANSFER  FCN 

3  Check  AK  input  selector  If  nit  2  channels 
inf if n  operator 

paisa  for  operator  -o  correct  it  ana 
loop  to  check  aqnn 
Clear  blocks  needed  for  averaging 


3  Zoon,  displaying  Channel  5  power  spectren 

2  Aver ace  spectra 

9  Loop  tor  specified  4  ;f  averages 

2  Conpote  transfer  function 

Sobrovtine  call  to  move  cata 
Retorn  to  calling  program 


SUBROUTINE  TO  N3Vt  DATA  BLOCKS 
Hove  data  blocks  fron  :oon  area 
to  area  for  operator  viewing 


Take  logs  of  spectra 


Retorn  to  calling  prsgran 

SUBROUTINE  FfliTofF  LINE”  ZOON  TRANSFi5^CN 
Clear  slocks  needed  -or  averaging 


2  Zoon  on  Channel  3 

1  Zoon  on  Channel  A 

Display  Channel  3  power  spectrvn 

2  Average  spectra 

9  Loop  tor  specified  *  of  averages 
2  Conpote  transfer  function 

Sobrovtine  call  to  »ove  data 
Retorn  to  calling  progrcn 

SBSSS*Sa3SS*3*SSSS3SS«tf28SSS3«SSSSSSSS*2SXSS 
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L 

-S7 

STACK  S’7 

1113 

POUEO  jPECTRM  PROGRAM 

Us 

37 

Recall  preuiius  ucriacle  parameters 

ft s 

17 

Y 

3828 

51 

Read  tot  Suffer  ISl 

v  U 

4 

1 

Print  nessajes  4  4  5 

y  a 

S 

l 

as 

512 

Bloc*  s»:e  512— user  *»«■•  cnange  after  pause 

r 

Pause  fur  operator  ac-ion 

Y  li 

6 

i 

Y 

i 

0 

Parameter  entry 

1  ■ 

2 

1 

Rifles;  *  •  1 

J 

110 

0 

Variable:  :4 

Y  . 

14 

29109 

1100 

BLOCK  1103 

Y  A+ 

0 

illl 

14D 

J 

19 

L 

lid 

4 

li-J? 

u 

1102 

»/ 

:0C 

0 

-2  Tec  2  *il*ers—auto  select  ncae 

L 

110? 

END  OF  Block  hog 

Y  U 

7 

1 

9 

Pause  for  operator  ic'icn 

Y 

5838 

52 

Y  J 

Id 

1 

RA 

c 

Analog  input  ia  'REPEAT'  tu  ror.it or  inputs 

r 

TL 

Pouer  ipect^u.'t  •  log >  for  user  reference 

; 

2100 

1 

1  Junp  t«  next  stack  ana  continue 

143 

i 
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STACK  SB 

LOOP  POINT  FOR  POWER  SPECTRUM 


L 

L 

Y 
1 

Y 

Y 
J 
r 
r 

Y 

Y 
J 

Y 
L 


A* 


5- 


-58 

2009 

5838 

11 

1 

2 

111 

15 
12 

1 

2 

111 

16 

1251 

1 

ID 

1251 

22 

3100 

1259 

1252 

13 

1 

2 

100 

17 

5838 

15 
1 
2 

101 

18 

1300 

0 

ID 

1301 

16 
1 
2 

100 

19 
17 

1 

2 

100 

20 

1309 

1302 

1303 


S2 

1 

1 

32767 

0 

211  ID 

1 

1 

2 

0 

2090D 
1251  16D 


1 


1 

1 

3 

0 

2000D 

S3 

1 

1 

2 

0 

2000D 
1300  18D 


i 

32767 

0 

2000D 

1 

1 

19D 

0 

2800D 

1 

1 


Paraneter  entry 
Range:  i  -  32767 
Variable;  IS 


Parameter  entry 
Range:  1  -  2 
Variable:  16 

BLOCK  1250  (ends  in  stack  60) 
Branch  oer  variable  16 
te  L12S1  or  1252 


Pause  for  operator  action 
Svbrovtine  call  -  baseband  neasurenent 
then  go  to  end  of  block  1250 


Parameter  entry 
Range;  1-3 
Variable:  17 


Paraneter  entry 
Range:  1  -  2 
Variable:  18 

BLOCK  1300  (ends  in  stack  59) 
Branch  per  variable  18 
to  L1301  or  1302 


Paraneter  entry 
Range:  1  -  32767 
Variable:  19 


Paraneter  entry 

Range:  1  -  (valve  of  var  19  -  1) 
Variable:  20 


STACK  S? 


■  I’1'  I  11  I  'J 


L 

-59 

1 

L 

1383 

y  u 

18 

1 

/. 

2888 

Y 

2883 

2881D 

Y  8 

19 

1 

/. 

2888 

/. 

-1 

Y 

2884 

28818 

L  V 

1358 

Y  IF 

2803 

28040 

2  -2 

Y 

2884 

20038 

Y  _ 

2883 

20018 

L  * 

1359 

Y  A- 

2888 

28848 

20838 

Y 

28 

28088 

Y  T 

2808 

Y  A* 

2801 

20038 

20888 

Y 

19 

28818 

L  ‘ 

1488 

Y  IF 

28 

198 

1  -2 

Y  A- 

28 

198 

1 

L 

1489 

L 

1389 

L 

1458 

Y  A* 

6 

1458 

178 

J 

80 

L 

1451 

Y 

5838 

54 

Y  V 

22 

1 

D 

J 

1459 

1 

1 

L 

1452 

L‘ 

1453 

Y  V 

20 

1 

Y 

1 

1 

Y  " 

2 

2 

; 

108 

0 

Y 

24 

28888 

L  " 

1508 

Y 

5838 

S4 

Y  A* 

8 

1588 

24D 

J 

9D 

L 

1581 

» 

V 

1509 

1 

t 

L 

1582 

Y  U 

28 

i 

Y 

l 

i 

Y  * 

2 

3 

T  “ 

0* 

100 

8 

L 

1556 

Y  A* 

8 

1558 

'28808 

J 

6D 

L 

1551 

D 

J 

1551 

L 

1552 

J 

L 

ffi 

1 

1 

J 

1555 

1 

1 

Cursor  on,  Cursor  parameters  to 
rariables  2800  -  2082  (2001=frequency) 
Variable  2803  gets  1st  cursor  frequency 

Cursor  on.  paraneters  to  2880  -  2682 
Corser  off 

Variable  2084  gets  2nd  cursor  frequency 
BLOCK  1350 

If  1st  cursor  frequency  >  2nd 
then  swap  then 
END  Of  BLOCK  13S8 


Variable  28  gets  :oin  bandwidth 


Variable  19  gets  zoon  center  frequency 
BLOCK  1408 

If  bandwidth  not  (  ctr  freq 
then  bandwidth  =  ctr  freq  -  i 
END  OF  BLOCK  1400 
ENB  Of  3L0CK  1380  (fran  stack  53) 

BLOCK  1450  (ends  in  stack  60) 

Branch  per  variable  17 
tu  L1451,  1452  or  1453 


Pause  for  operator  action 


Paraneter  entry 
Range:  1-2 
Variable:  24 

BLOCX  1506  (ends  in  stack  60) 

Branch  per  variable  24 
tu  L1541  or  1502 


Paraneter  entry 
Range:  l  -  J 

Variable:  2080  itenporarv) 
BLOCX  1558  (enas  in  stack  60) 
Branch  per  variable  2008 
to  LiSSi,  1552  or  1553 

Pause  6  loop  to  here—abort 
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w  -- . — .  -  — - ^ .a-,  r  — dMflMiiifci  mii  ■  n  r  i 


STACK  60 


l 

l 


-60 

1S54 


Y 

22 

134 

Variable  22  gets  default  start  track 

Y  ' 

J  " 

21 

1559 

64 

variable  21  gets  default  0  of  records 

L 

1555 

y  y 

27 

1 

?  - 

i 

19^ 

Paraneter  entrv 

Range:  1  -  197 

J  ■ 

100 

0 

Variable:  22 

Y 

22 

20000 

Variable  22  gets  user's  start  track 

Y  «- 

21 

198 

20000 

Variable  21  gets  0  of  records  left 

END  OF  BLOCK  1550  (fron  stack  59) 

L 

1559 

MS 

32 

220 

Position  throughput  file  to  start  track 

Y  y 

22 

i 

D 

Pause  for  operator  action 

Y  BS 

13 

Variable  l3  aets  current  block  size 

BS 

4096 

Set  block  size  to  nax  for  throughput 

MS 

22 

1 

210 

Perforn  ADC  throughput 

BS 

i3D 

Restore  block  size 

y  y 

21 

1 

L 

1509 

END  OF  BLOCK  1500  (fron  stack  59) 

L 

1459 

END  OF  BLOCK  1450  (fron  stack  59) 

L 

1600 

BLOCK  1600 

Y  A* 

8 

1600 

170 

Branch  per  variable  17 
to  LloOl,  1602  or  1603 

J 

OS 

L 

1601 

Y 

41 

190  200 

Initialize  zoon— on  line  preprocessor 

J 

4100 

1 

i 

Subroatine  call— zoon  neasurenent 

Y 

40 

0 

Reset  zoon  to  baseband 

J 

1609 

Ge  to  end  of  block  1601 

L 

1602 

Y 

43 

190  200 

220 

Initialize  zoon— off  line  preprocessor 

J 

4500 

1 

1 

Subroutine  call— zoon  neasurenent 

Y 

43 

0 

Reset  zoon  to  baseband 

J 

L 

T 

I 

Y 
L 
L 

Y 

Y  U 

Y  U 

Y  U 
MS 
MS 
D 

J 


1619 

1683 

41  19D  200  22D 


4S0I 

41 

1609 

1259 

5838 


37 

27 

2101 


23 

25 

26 


1 


54 


1 


-2 


Go  t«  end  of  block  1601 

Initialize  zoon— off  line  software 
Subroutine  call— zoon  neasurenent 
Reset  toon  to  baseband 
END  OF  BLOCK  1600 
END  OF  BLOCK  1250  (fron  stack  58) 

Power  spectrvn  conplete 
Print  nessages  to  operator 

Save  variable  paraneters 

Pause  for  operator  action 
Ljop  to  repeat  power  spectrun 
-  - SPEgTRUft  (tgCfpt  subroutines) 
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STACK  si 

BASEBAND  POWER  SPECTRUM  SUBROUTINE 


L 

-61 

L 

3890 

Y 

5833 

55 

Y 

6 

i 

3 

Y  W 

32 

1 

D 

J 

3888 

-1 

CL 

i 

L 

3810 

RA 

9 

1 

Hi 

8 

F 

a 

0 

8 

SP 

9 

t 

3910 

15D 

0 

j 

9055 

1 

1 

t 

i 

16 

; 

i 

3 

x< 

t 

TL 

9 

< 

L 

4890 

Y 

5838 

55 

Y 

6 

1 

3 

Y  i 

32 

k 

D 

J 

4098 

-1 

CL 

2 

L 

4810 

Y 

45 

i 

l 

SP 

1 

* 

4810 

1SD 

9 

x< 

2 

TL 

<  . 

4500 

1  OU1 
.•h-t 
in 

i 

i 

D 

2 

4 

SP 

i 

* 

4510 

1SD 

8 

X< 

2 

TL 

Check  ADC  input  selector.  If  not  1  cnannei 
infer*  operator 

pause  for  operator  -o  correct  it  and 
loop  to  check  again 
Clear  block  needed  for  averaging 

Analog  input,  displaying  average 
Hann 

Transfer* 

Clear  DC 
Average  spectrun 

Loop  tor  specified  *  of  averages 
Subroutine  call  -  clear  2nd  half  block 
Hanning  correction  for  broadband  noise 

Load  average  to  block  3 
Take  log 

Return  to  calling  progran 

SUBROufl«~FO»-ONLINE~foidhTGkeER~SPEcfRUH 

Check  ADC  selector.  If  not  2  channels 
infer*  operator 

pause  for  operator  to  correct  it  and 
„  loop  to  check  attain 
Clear  block  needed  for  averaging 

Zoon,  displaying  pouer  spectrun 
Average  spectrun 

Leap  tor  specified  ♦  «f  averages 
Load  average  to  block  I) 

Take  log 

Return_to  calling  progran 

SUBROUTINE  FOR  OFF  LINE  ZOOM  POWER  SPECTRUM 
Clear  block  needed  for  averaging 

Zoon 

Displav  average  (i  sweep) 

Average 

Loop  for  specified  f  of  averages 
Load  average  to  block  3 
Take  log 

Return  to  calling  progran 


L 

-62 

L 

9151 

9057 

Y 

29 

; 

9056 

-l 

L 

9157 

25D 

a 

4 

a 

5 

2SD 

i 

9159 

-V 

L 

9155 

r 

29 

9058 

j  ’ 

9056 

-i 

L 

9058 

2SD 

a 

i 

j 

9059 

-l 

L 

9056 

r  a* 

0 

9951 

j 

0D 

-i 

L 

9052 

r  bs 

i3 

Y  : 

26 

13D 

Y  ; 

25 

260 

J 

29D 

-i 

L 

9051 

L 

9059 

< 

CLEAR  2ND  HALF  BLOCK  SUBROUTINE 
Set  retirn  label  valve 

Retvrn  label 

26D  Clear  last  half  black  4 
260  Clear  last  half  block  5 

Entry  pt.  Pur  Spec  pqn 
Set  return  label  valve 

Retirn  label 

26D  Clear  last  half  block  1 


14D  Has  HP  filter  selected? 
No,  Gto  9951 

Yes 

Get  current  blocksne 
Store  Bs/2  in  variable  2b 
Store  Bs/4  *n  variable  25 
Indirect  retvrn 

End  of  pqn  905049055 
Retirn  to  calling  prograft 


Test  buffer  nessages  for  Transfer  Function  and  Pouer  Spectrun  prograns 
BUFFER  0  MESSAGE  * 

Si  1 

HP  TRANSFER  FUNCTION  PROGRAM 


51  2 

SELECT  EXCITATIGN  TYPE 
51=RAND0M  -  BASEBAND  OR  ZOOM 
2=HP  DAC  -  BASEBAND  ONLY 
3=TRANSIENT  -  BASEBAND  ONLY 


Si  3 

INPUT  DESIRED  DAC  OUTPUT  IN  W 


Si  4 

HP  POUER  SPECTRUM  PROGRAM 


ST  ADC  FREQUENCY  RANGE  AS  DESIREI 
(SAMPLE  MODE  4  MULTIPLIER) 

SET  ADC  TRIGGER  TO  'FREE  RUN 
CHANGE  BLOCK  SIZE  IF  DESIRED 
PRESS  ‘CONTINUE’  WHEN  READY 


51  6 

ARE  HP  FILTERS  INSTALLED? 
Q*NO,  Si=YES 
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51  7 

SET  KEYBOARD  REPEAT/SINGLE  SWITCH 
TO  'REPEAT' 

PRESS  'CONTINUE'  WHEN  READY 

51  98 

UPPER  UNIT  = 

51  99 

ENTRY  OUT  OF  RANGE-PLEASE  REENTER 
LOWER  LIMIT  = 

PRESS  CONTINUE  WHEN  READY 


52  8 

TURN  ON  RANDOM  EXCITATION  SOURCE 


52  9 

SET  TRIGGER  SOURCE  AS  DESIRED 
WPACT  STRUCTURE  REPEATEDLY 


52  1# 

SET  OVERLOAD  VOLTAGES  AND  TRIGGER 
LEVELS  FOR  SIGNAL  AMPLITUDES 
NOVc  fEPEAT/SINGLE  SWITCH  TO 
'SINGLE*  WHEN  READY.  IF  SOURCE 
IS  NOT  IN  FREE  RUN.  TRIGGER  THE 
SYSTEM  AGAIN  TQ  CONTINUE. 


S2  ii 

ENTER  NUMBER  OF  AVERAGES  DESIRED 


ENTER  hIaSUREMENT  TvPE 
’“BASEBAND 
S2=ZQ0M 


52  13 

ENTER  ZOOM  MEASUREMENT  MODE 
i=QN  LINE.  PREPROCESSOR 
S2=QFF  life,  PREPROCESSOR 
3=OFF  LINE,  SOFTWARE 
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52  14 

ZOOM  NOT  APPROPRIATE  WITH  HP  OAC 
BASEBAND  MEASUREMENT  WILL  BE  HADE 


53  15 

HOW  HILL  YOU  SPECIFY  ZQOH  BANDWIDTH? 
i=NUM£RIC  ENTRY  -  CTR  FREQ  i  BU 
2*CURS0R  -  ON  PRIOR  NEASUREHENT 


S3  16 

ENTER  CENTER  FREQUENCY 


53  17 

ENTER  BANDWIDTH 


53  18 

HOVE  CURSOR  TO  START  FREQUENCY 
PRESS  'VALUE'  (SWITCH  REGISTER  11) 


53  1? 

HOVE  CURSOR  TO  END  FREQUENCY 
PRESS  'VALUE* 


53  20 

ANALYZE  OLD  OR  NEW  DATA? 
i=QLD  (FROH  THROUGHPUT  FILE) 
2=NEW 


54  21 

THROUGHPUT  COMPLETED 


54  22 

PRESS  'CONTINUE*  FOR  MEASUREMENT 


54  23 

MEASUREMENT  COMPLETE 


54  24 

TO  DISPLAY  RESULTS.  PRESS: 
'DISPLAY'  '0*  LliG  TRANSFER  FCN 
•DISPLAY'  *1'  COHERENCE 
■DISPLAY*  '2*  INPUT  POWER  SPECT 


I 


•DISPLAY* 


OUTPUT  POWER  5PECT 


•DISPLAY*  *S4*  CROSS  POWER  SPECT 


54  25 

TO  COPY  DISPLAY  ON  TERMINAL: 
PUT  TERMINAL  IN  GRAPHICS  MODE 
PRESS  *GOU)  KEY*  ‘PLOT* 


54 

TO  MAKE  ANOTHER'/iEASJRSMENT: 
PUT  TERMINAL  IN  ASCII  MODE 
PRESS  ‘CONTINUE* 


ENTER  STARTING  TRACK  FOR  THROUGHPUT 


54  20 

THROUGHPUT  WILL  USE  TRACKS 
135  THROUGH  198  ON  THE  LOWER 
(FSDS)  DISC.  IS  THIS  OK? 
i=NQ  -  ABORT 
2= YES  -  PROCEED 
3=NQ  -  ASK  ME  FOR  TRACK  ♦ 


55  29 

IMPACT  STRUCTURE  ON  CUE  (BEEP) 
FOR  EACH  AVERAGE 


55  31 

SET  ADC  INPUT  SELECTOR  TO  "A8* 
PRESS  'CONTINUE*  WHEN  READY 


Buss  rs  cot  operator 
Beeps  to  :«  operator 


55  32 

SET  ADC  INPUT  SELECTOR  TO  *A* 
PRESS  'CONTINUE*  WHEN  READY 


Beeps  to  cue  operator 


APPENDIX  C 


CALIBRATION  DATA 


ACC&EROMETER  MOUNTING  WAX 
Model  080A24 


•W20TH0NICS 


Petro-Wax  functions  to  transfer  motion  from  the  test  object  to  the  sensor.  It  is 
used  to  couple  the  sensor  directly  to  the  test  object.  It  is  a  convenient,  temporary  mount¬ 
ing  replacement  for  studs  and  permanent  adhesives. 

An  elastic  interface  structure  with  ticavy  damping,  it  forms  quickly  to  irregular  sur¬ 
faces  to  facilitate  easy  mounting  of  sensors.  This  pliable  uax  enables  the  sensor  to  be 
mounted  in  nearly  any  convenient  spacial  coordinate  system. 

In  conjunction  with  the  sensor,  the  vox  forms  a  heavily  damped  spring-mass  system. 

The  frequency  response  is  a  function  of  transducer  mass,  mounting  area,  depth  of  wax,  and 
test  temperature.  As  the  amount  of  uax  at  the  interface  of  the  sensor  and  test  object  in¬ 
creases,  the  first  resonance  of  the  system  decreases,  which  limits  the  frequency  response 
of  the  fixturing  (see  graphs).  This  mrunting  technique  is  primarily  for  use  at  room  temp¬ 
erature.  The  wax  cannot  be  used  effectively  at  high  or  low  temperatures. 

How  to  use:  (1)  To  insure  a  secure  and  valid  fixture,  be  certain  all  surfaces  are  free 
ot  oil  and  dirt. 

(2)  Apply  the  wax  directly  to  the  base  of  the  transducer  or  to  the  adhesive 
mounting  base.  The  amount  of  wax  used  depends  on  your  individual  appli¬ 
cation.  Use  of  the  adhesive  mounting  base  helps  keep  the  transducer  clean. 

(3)  Press  the  sensor  firmly  against  the  test  structure  to  insure  secure  mount¬ 
ing  and  as  little  uax  at  interface  as  possible.  This  provides  the  best 
frequency  response. 

(4)  Precede  with  measurement. 

PETRO-WAX  is  available  in  quantity  from:  Katt  &  Associates 

P.0.  Box  98269 
Pittsburgh,  PA  15227 
(412)  885-5727 


MODEL  NO.  308° 
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MODEL  NO 


30ZA 


Ihin  <0.010' 


Horizontal  Scale  10  000  Hz 
Mag.  :  wert.  scale  10  db/diw 
Phase:  wert.  scale  50  a/diw 


Thin  Med  Thick 


fREO 

(Hi) 

oev 

(S> 

ocv 

W 

DCV 

(%) 

20 

-l 

-l 

-i 

SO 

0 

0 

-i 

100 

0 

0 

0 

200 

0 

0 

0 

500 

1 

1 

2 

1  000 

1 

2 

8 

2  000 

1 

T 

41 

S  000 

1 

62 

35 

10  000 

11 

56 

-450 

20  000 

Magnitude 

_ _ 

! 

Phase 

f  rcqucncy  — ► 

_ L.  1  -  -  t  -  1  t  I 
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MODEL  NO.  30 J A 

- 1 - 1 - 1 - 1 - 

Thin  00.010" 

Maqnitude 

T - 1 - 1 - 

i - 1 - 

Horizontal  Scale  25  000  Hz 

~ 

1 

Mag.  :  vert,  scale  10  db/div 

Phase 

Phase:  wert.  scale  50  °/diw 

. .  J _ ! _ 1 _ 1 _ 

!  1 

(  25  KHz  ^ 

Thin  Med  Thick 


FREQ 

oew 

DCV 

DEV 
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(S) 

<■;> 
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20 

50 
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100 
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200 
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0 

500 
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11 

A8 

10  000 
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20  000 
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rlB  CALIBRATION  CERTIFICATE 
IMPACT  FORCE  HAMMER 

Model  Ho.  0863  03 
Transducer  Model  No.  — 

Transducer  Serial  No.  9  _ 

Hammer  Calibrator:  PCB  ModeJ_ 
Pendulous  Mass  0-79- lb  [jpCr  gm) 


I  LEU  FIEZCTRCNICS.  INC 
Customer: 

Invoice  No: 


Cate: _  ?~a2  3  ~ 


including  accelerometer 


Initials: 


Basic  Force  Transducer  Sensitivity  mV/lb  (reference) 

(for  stationary  installations) 

Reference  Accelerometer  Sensitivity  /£/£-  mV/q  (302At*^/N  75~3/) 

Ratio:  Force  Transducer  Stnsi tivi ty/Accelerometer  Sensitivity_  /-or,0. 

(from  test  of  mass  impacting  stationary  banner) 

HAMMER  SENSITIVITY  (3)  (with  2.7  gm  Al/Plastic.  Al/Steel,  Al/Rublor  T:ps 


Configuration 


with  with  with 

Steel  Extender  Al.  Extender  No  Extender 


Ratio:  Force/Acceleration  /  / 

Sensitivity  (2)  /.  uj- 

Hammer  Sensitivity  nV/lb 


Difference:  (1) 


-J./Z 


’P^Z 


NOTES:  (1)  Difference  from  reference  sensitivity  of  force  transducer 

(2)  In  transfer  fuction  testing,  the  important  factor  is  the  ratio  of 
sensitivities. 

(3)  Because  of  normal  behavior  of  the  banner  structure,  the  aoparent 
sensitivity  of  the  haunter  in  motion  differs  from  the  stationary 
calibration  of  the  force  transducer.  It  is  less  by  a  factor 
proportional  to  the  ratio  of  the  "iass  of  the  impact  cap  and  seismic 
distributing  mass  in  the  transducer  to  the  total  mass  of  the 
hanmer  structure.  Using  a  heavier  harrier  head  or  installing  the 
force  transducer  on  the  structure  only  changes  the  problem.  A 
heavier  harrier  head  tends  to  penetrate  the  test  object  or  cause 
multiple  bouncing.  When  mounted  on  tne  test  object,  the  inertial 
mass  in  the  transducer  causes  it  to  act  as  an  accelerometer  sensing 
motion  of  the  test  object. 


t.  VOLTAGE  SENSITIVITY  . 

2.  MAXIMUM  TRANSVERSE  SENSITIVITY 

3.  RESONANT  FREQUENCY  . 

4.  DISCHARGE  TIME  CONSTANT  . 

5.  OUTPUT  BIAS  LEVEL  . 

6.  FREQUENCY  RESPONSE: 


I  Mm  Temp 


mv/g  @  1 00  Hi.  8  g  s  PK 


percent 

KH, 

— 


CNiferation  treceeble  io  NBS  through  project  no. . 

»CB  PIE20TTONICS  INC 
3425  WALDEN  AVENUE 
OEPEW  NEW  YORK  14043 


00  iV>. 3 


77  7/ 


3000  5000 


w-1 
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K3  CALIBRATION  CERTIFICATE 
IMPULSE  FORCE  HAJ^ES. 


5^  t  b  •  J  ^  i  (• 


ii  PIEZOTRCNICS,  IRC. 


fatal  Ns.  086PgQ 
Serial  He.  23S 


KH  Ho *€/rj?//7<g* 


CklOtOMtl 

A&sat 


hm  lb 

Linearity  arret  ■*-  e2.°  % 

Discharge  Tim  Constant  ^  a 

Output  laaedance  /c’*'  ohas 

Output  Blaa  /<?•  J.  valta 

Iracaabla  to  NBS  throu^i  7S~? J 

Initial*.  Qatar  V  / £2 


Invoice  No.i 


Aoaa  itron*  ten  Nodal  No.  iCf&S  .  Serial  Wo.  Sana.  2  «V/a 

Pendulous  Toot  Haas  —  lb  (  c.raa)  including  accalaroaater 

Sensitivity! 


COWICURATIW 

-  '  " 

Tip 

STiCL 

Vi  M>fL 

Catandar 

Ncalt 

S76CL 

SCALING  r ACT OR 
(SENSITIVITY  RATIO) 

,ns> 

.  \z4 

HMfCft 

■V/lb 

U3.3 

7o.  9 

SENSITIVITY 

(aWAO 

zss 

/S.9 

WfCSt 

1.  The  aanaitivity  ratio  (Sa/Sf )  ia  the  seal log  factor  for  converting  structural  transfer 
aesauremnta  into  anglnsaring  units.  Divide  results  by  this  ratio. 


2.  Each  specific  hooter  configuration  has  a  different  sensitivity.  The  difference  is 
a  constant  percentage,  ohicti  depends  m  the  asaa  of  the  cap  and  tip  asianbly  relative 
to  the  total  aaca  of  the  head.  Calibrating  the  specific  heaiar  structure  being  used  euto- 
aatlcally  canpenaates  far  ease  effects. 


Effective 


ossa  £ 


with  X2A07  attached  and  vinyl  capped  plastic  tip. 


■I  > 
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ulftOjgr  ,lt* 


i 


■aHUMMAMI 


mnwadmit  Pr,ib(^r.  ^  -r*. 


tor  No.  NllZSmZltA  Hid. 


HD 


rrp 


uu 

PtCZOTROMICS 


CALIBRATION  DATA 

(o> 

I.  C.  P.  ACCELEROMETER 
(por  I S  A  S37J2) 


fl  I 'o  -t-WwioAu^" 


Medal  Na.. 

Send  Ns.  _ 

Rmtt 

Mailapat. 

/*** 

|*S 

MasTseie 

'  s'** 

*F 

1.  VOLTAGE  SENSITIVITY 


£6L 


mv/g 


2.  MAXIMUM  TRANSVERSE  SENSITIVITY 

3.  RESONANT  FREQUENCY  . 

4.  DISCHARGE  TIME  CONSTANT  . 

6.  OUTPUT  BIAS  LEVEL  . 

6.  FREQUENCY  RESPONSE: 


<P.8 


<9.X 


12_ 


9  100  Hz.  8  g's  PK 

percent 

KHe 

seconds 

Volts 


Frtg.  Hz 

■ 

1  10  i 

30 

so 

100 

300 

500 

1000 

3000  { 

5000  | 

7^ 

/00CO 

Deviation  % 

EH 

EH 
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Calibration  traceable  to  NBS  through  project  no.  _ 
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1.  VOLTAGE  SENSITIVITY 
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2 .  MAXIMUM  TRANSVERSE  SENSITIVITY 

3.  RESONANT  FREQUENCY . 

4.  DISCHARGE  TIME  CONSTANT  . 

6.  OUTPUT  BIAS  LEVEL  . 

8.  FREQUENCY  RESPONSE: 
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Calibration  traceable  to  NBS  through  project  no. . 
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I.C P.  TRANSDUCER  DATA 
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Sensitivity _ ^ , 

Linearity  *-/•* 
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Cal.  Range  - BUFFALO,  HEW  YORK  H22C 

Input  Time  Constant  ■>***  >«t  y 
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LOCALLY  GENERATED  USER  KEYBOARD  PROGRAMS 

The  locally  generated  keyboard  programs  presented  in  the 
next  six  pages  allow  the  use  of  impulse  technique  or  random 
excitation  while  utilizing  the  Modal  package  of  the  HP-5451C 
Fourier  Analyzer.  Each  program  and  subprogram  is  identified 
by  the  first  line  of  the  program  listing.  Program  number  1 
(identified  by  the  -1  on  line  1)  is  the  control  program, 
and  from  this  program  subprograms  50,  51,  52,  53,  54,  58 
and  59  are  called  as  needed. 

All  eight  programs  are  automatically  loaded  into  the 
computer  memory  when  the  following  command  is  executed,  as 
discussed  in  Chapter  V,  step  number  35, 

[JUMP]  0  1  [ENTER] 
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TEST  CHAMBER  CHARACTERIZATION  DATA 
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SETUP! TEST  I DATE)  TYPE  1  RANGE  I INPUT  LOC.I  RESPONSE  1 FILTER t STORE  I  DATA 
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SETUP t TEST  I DATE  I  TYPE  I  RANGE  I INPUT  LOC.  I  RESPONSE  I FILTER t STORE  I  DATA  I 

NO.  I  NO. I  | (BorZ) | (MAX. FREQ) I (GRID  NO.) I  LOC.  |(YorN>|  DATA  I  LOC.  I 
III  I  or  I  ((GRID  NO.) I  I (Y .N) I (RECORD) I 


43PTT 


Jk066 

tfoTl 

P03> 

g  063 


I  }<!/  i  JCPH/O*.  J 

*L_  fe. T !  Ao  _  £  &£JL9so !_!- %  6 

1 

L09\ 

i  .  i  .  i 

IJL.IJL! 

2 

.go?_L 

DoTt 

6063 


ICF  %Z<\  I 

Icf-Ck.tl1.!__upja6 


_<o£L 

_j>o_S>„ 

6063 


UftQ  Ftcnr  ReMo^ep 

T&*fT  :  3*i°  C 


169 


fit 


ve^pr:  2^°  c 


CH*WKT£V?  I  -frA-TI  C** 


AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 
CNTR  FREQ  llAl-  H*  DELTA  FREQ  ^Sl^Hx 

set  up  no.  2L3« _ i.l.  40A*_  p.l.  CASlL- 


f7- 


NAT.  FREQ 
MODE  (Hz) 


DAMP.  FACT. 
(X) 


DAMP.  COEFF. 
(RAD/SEC > 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 

CNTR  FREQ  2*111  Hz  DELTA  FREQ  £\£*£hz 
SET  UP  NO.  gjtf  __  I.L.  4 Qhk-  P.L.  h-QSil L 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 

CNTR  FREQ  Hz  DELTA  FREQ 

SET  UP  NO.  XLS___  I.L.  43$>A_  P.L. 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 

CNTR  FREQ  35 Hx  DELTA  FREQ^*A17hz 
SET  UP  NO.  _ I.L.  P.L.  <*S>SZL 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 

CNTR  FREQ  *jtQ2s.  Hz  DELTA  FREQ  /j.9S?  Hz 
SET  UP  NO.  S_"7 _ I*L.COfc6 _ P.L.  4 PX1_ 
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AVERAGE  MODAL  FREQUENCIES  ANO  DAMPING 

CNTR  FREQ  Hz  DELTA  FREQ  AJWhz 

SET  UP  NO.  I.L.  P.L.  LOT/ 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 


CNTR  FREQ  SMi 
SET  UP  NO.  225. 


H*  DELTA  FREQ  44?Ahz 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 


CNTR  FREQ 
SET  UP  NO. 


>_7  Hz  DELTA  FREQ  Qj__ 
LQ—  I.L.  P.L.  4A. 
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DAMP.  FACT.  DAMP.  COEFF. 
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PRELIMINARY  MODE  IDENTIFICATION 

CNTR  FREQ  7 f*~  Hx,  DELTA  FREQ  t'A2?-  Hz 
SET  UP  NO.  _£_£_•  I.L.  P.L.Ji.3iC£ 

NODES  INTERVAL  LEVEL  ITER. 

-2  O  /O.O  i 

NODE  FREQ.  OAHP  AMPL.  PHS 
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PRELIMINARY  MODE  IDENTIFICATION 
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PRELIMINARY  MODE  IDENTIFICATION 
CNTR  FRED  2idJr  _  Hz.  DELTA  FREQ  Hz 
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PRELIMINARY  MODE  IDENTIFICATION 
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PRELIMINARY  MODE  IDENTIFICATION 
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PRELIMINARY  MODE  IDENTIFICATION 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMRING 
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SET  UP  NO.  Jlj _ I.L.  _ P.L.  _ 
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AVERAGE  MODAL  FREQUENCIES  ANO  DAMPING 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 
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AVERAGE  MODAL  FREQUENCIES  AND  DAMPING 
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PRELIMINARY  MODE  IDENTIFICATION 
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